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COMPLETE  REPORT 
Summary 

The  non- destructive  control  systems  with  electromagnetic  UHF  probe  of  plane-layered  structures  are 
mainly  used  for  the  estimation  of  the  location  of  defect.  As  a  rule  they  work  by  a  given  frequency  and  can¬ 
not  determine  neither  the  depth  of  the  defect  location,  nor  its  proportions.  Thus,  these  are  only  defecto- 
scopes  aimed  at  locating  defects.  There  also  exist  systems  oriented  at  the  estimation  of  plane-layered  di¬ 
electric  structures’  parameters,  but  within  a  very  narrow  diapason  of  their  change.  The  existing  systems  do 
not  ensure  the  control  and  diagnostics  (defectoscopy  and  defectometria)  of  multi-layered  constructions  in 
the  wide  diapason  of  the  change  of  their  geometrical  and  electro-physical  parameters. 

In  this  report  we  propose  a  technique,  which  allows  estimating  the  amplitude  and  the  phase  of  the 
complex  reflection  coefficient  for  various  frequencies.  At  the  beginning  stage  of  the  location  of  defect  the 
device  operates  by  a  given  frequency,  similar  to  defectoscope.  Under  such  working  conditions  the  zone  on 
the  surface  of  the  plain-layered  structure,  under  which  the  defect  is  possibly  located,  is  found.  Next,  the 
mode  of  multi-frequency  sounding  of  the  given  area  comes  into  operation.  That  is,  in  the  area  where  the 
defect  is  located  the  device  estimates  the  complex  reflection  coefficient  by  various  frequencies  (these  coef¬ 
ficients  bear  the  information  about  the  dielectric  structure  with  a  defect.)  The  information  about  the  value 
of  complex  reflection  coefficient  is  saved  in  the  memory  of  the  computer.  After  that,  the  mathematical 
technology  comes  into  operation.  It  serves  for  the  estimation  of  the  parameters  of  the  defect  (the  depth  of 
the  location  of  the  defect  and  its  thickness).  The  model  of  interaction  of  the  electromagnetic  wave  with  the 
multi-layered  structure  plays  the  important  role  in  this  approach.  Using  such  model  the  complex  coeffi¬ 
cients  of  reflection  by  the  necessary  frequencies  are  estimated  by  means  of  the  mathematical  procedures. 
Next,  the  functional  are  built,  which  , compare  experimental  and  theoretical  complex  reflection  coefficients 
in  such  a  way  that  the  values  of  parameters,  that  ensure  reaching  of  the  global  minimum  of  this  functional, 
are  equal  to  the  values  of  the  parameters  of  the  defect  or  of  the  parameters  of  the  layers  in  the  structure. 

Here  we  give  the  complete  report,  which  includes  the  methodology  of  resumption  of  parameters  of 
plane-layered  dielectric  structure  basing  on  the  comparison  of  the  theoretical  and  experimental  data.  The 
designing  of  the  radio-wave  NDT  device  for  dielectric  structures  has  finished.  Frequency  range,  sensitiv¬ 
ity,  measurement  precision,  models  of  measuring  system  and  interaction  of  probing  field  with  dielectric 
flat-layered  environment,  as  well  as  methods  for  solution  of  inverse  problems  would  allow  solving  a  wide 
range  of  practical  problems  of  non-destructive  testing. 
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I.  DEVELOPMENT  OF  THE  MATHEMATICAL  MODEL  FOR  FIELD  CALCULATION  IN  THE 
ANTENNA’S  APERTURE,  WHICH  TAKES  THE  ANTENNA’S  GEOMETRY  INTO  ACCOUNT 


1.  Introduction 

Contemporary  high-quality  analysis  of  the  field  of  horn  antennae  in  the  radiation  zone  is  based  on 
the  usage  of  the  solutions  of  problems  of  electromagnetic  excitation  of  semi-infinite  perfectly  conducting 
cones.  The  solutions  of  the  wave  diffraction  problem  by  the  semi  infinite  cone  are  well  known  (see  [1,2]). 
However,  they  do  not  take  the  diffraction  effects  on  the  edges  of  finite  cone  into  account.  These  effects 
lead  to  the  complication  of  field  in  the  antenna’s  aperture  and  losses  of  energy  because  of  the  penetration 
of  the  field  into  the  shadow  zone.  As  the  near  field  is  important  in  our  case,  the  consideration  of  these  ef¬ 
fects  or  their  minimization  is  one  of  the  key  problems.  The  aim  of  the  project  at  this  stage  is  to  show  the 
possibility  of  minimization  of  edge  effects  by  choosing  the  horn  antenna  parameters  (radius  of  aperture, 
length  of  the  cone’s  generate,  source  of  excitation).  The  real  field  of  horn  antenna  has  of  a  complex  charac¬ 
ter  and  is  formed  as  a  result  of  the  modes’  distribution,  their  transformation  and  interaction  on  the  edges, 
caused  by  the  vector  character  of  electromagnetic  field.  To  reach  the  goal  at  this  stage  of  the  project  we 
considered  the  simplest  case,  in  particular  the  formation  of  axially  symmetric  field.  Thus,  to  reach  the  aim 
we  will:  firstly,  solve  the  simplest  axially  symmetric  problem  of  radiation  field  of  the  radial  electric  dipole 
from  the  finite  conical  horn  antenna;  secondly,  we  estimate  the  field  of  the  adequate  semi-infinite  cone; 
thirdly,  for  every  scatterer,  for  various  combinations  of  their  geometrical  parameters  and  dipole  locations, 
we  find  numerical  distribution  of  the  fields  and  their  mode  structure.  The  data  obtained  will  be  used  for  the 
formation  of  the  principles  for  choosing  of  the  parameters  of  horn  antennae  and  the  method  of  their  excita¬ 
tion. 


2,  The  field  calculation  technique 

Let  us  consider  the  perfectly  conducting  conical  horn  antenna  of  finite  length  as 

Q :  {r  e[0,c];  0  <  6<  y;  <p& [0,2zr]} ,  (1) 


where  (r,0,cp)  is  the  spherical  coordinate.  The  horn  antenna  (1)  is  excited  by  the  radial  electric  dipole, 


which  is  located  at  the  conical  axis  in  the  point  with  the  coordinates  r  =  1  <c,  9  =  0  (see  Fig.  1).  In  this  case 
we  have  an  axially  symmetric  electromagnetic  field  with  nonzero  field  components  Er ,  E0 ,  H (p  (axially 
symmetric  TM-  waves).  The  time  dependence  of  the  electromagnetic  field  takes  as  exp (-icot)  and  is  omit¬ 


ted  in  this  paper. 

Using  the  semi-inversion  technique,  the  problem  of  the  elec¬ 
tromagnetic  waves  radiation  from  the  horn  antenna  ( 1 )  is  reduced  to 
the  solution  of  the  set  of  linear  algebraic  equations  of  the  second 


kind  as  in  [3,4] 

X  =  A~l{A- Ax)X  +  A~lF . 

(2) 

Here  X  =  {x„}”=1 

is  the  vector  of  unknown  coefficients;  A1 , 

A  are 

the  infinite  matrices  with  elements: 

P\W  K  E 

a(D_  L  -?JA 

(3  a) 

qn  (42-z„2)AZn(A)^(A) 

aqn  =(4  ~Z„yl 

(3b) 

with  q,n=  l,oo;  /VQ,  Av(0  are  the  modified  Bessel  and  Macdonald 


0=0 


Fig.  1.  Geometry  of  the  problem 
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functions  respectively;  p  =  -ike  k  is  the  wave  number,  c  is  the  horn  antenna’s  length;  zn  &{zp }"=1, 

e  \q  Yp=\  are  the  growing  sequences  of  the  positive  simple  zeros  and  poles  of  the  meromorphic  func¬ 
tion 


M(v) 


COS(JTV) 

(v2  -  0.25)PV_, , 2  (cos  Y )PV_! / 2  (- cos y) 


Pv_i/20  is  the  Legandre  function, 


W[KvIp ]p  =  Kv (p )/; (p)  -  K'v (P)  • 
The  matrix  elements  rnq  for  the  inverse  operator  A  1  are  defined  as  follows 

Tnq  =  (j[M_  (£q )]  1 1  [  M_ (z„ )] '  (z„  -f9)j  ■ 


(4) 

(5) 

(6) 


Here  [A/_(z)]'  -d / d\’{M  ( v)} v=_ ;  M.  (v)  is  the  split  function,  which  we  obtain  after  the  procedure  of 
factorization  of  the  function  (4).  The  known  vector  F  —  { is  given  as 

A  n  I^(P0)  (^-0.25)P4_1/2(-cosr) 


fk  =  - 


2\[f  LiSP i)  cos 0*4) 


(7) 


where  A0  =  2 P0Z //  ;  P0  is  the  moment  of  the  electric  dipole,  Z  is  the  wave  resistance  of  the  medium, 
P  0  =-ikl- 

Next  we  express  the  field  radiated  from  the  horn  antenna  through  the  solution  of  the  equation  (2)  as 
follows 


1  »  P  (cos  9)  K(p) 

U (r, 9)  =  -=  2>„  ,  .  - .  ,r>c,p=  -ikr  . 


fp  ”  (z2-0.25)P  (cosy)*^) 


(8) 


The  solution  of  the  similar  wave  diffraction  problem  for  semi  infinite  perfectly  conducting  cone  we 
represents  as 


U(r,0)  =  - 


A„n  ”  %^,„-i/2(cos(9)P„  _1/2(-cosy)  f^„(p)^„(P0)l  r>/, 


K  ”=1  cos(7rrjn) — rP_i/2(cosr)l  r<^’ 

L  '  Jv=nn 


(9) 


where  r\n  are  the  positive  roots  of  the  transcendental  equation  if  1/2(cos  y)  =  0 ;  0  <  0  <  y  . 

The  determination  of  the  electromagnetic  field  components  is  based  on  the  equations  (8),  (9)  and  differen¬ 
tial  operators  as 


[Er,Ee,H9) 


1  d  ,  .  ndU  1 

- (sin# - ); - 

r  sin  9  36  39  r  dr 39 


(, rU );  ikZ1 


(10) 


We  have  analysed  numerically  the  near  horn  antenna  fields  as  well  as  far  radiated  fields’  pattern.  The  far 
fields’  pattern  characteristics  were  calculated  by  the  formula 

D(9 )  =  lim|  rH(pe~ikr\  ■  (H) 

r— >oo 

Due  to  axial  symmetry  of  the  problem,  plots  of  diagrams  are  drawn  in  the  angle  range  0  <  9  <  n  . 

We  also  analysed  the  total  and  diffracted  fields.  The  expression  for  diffracted  field  potential  looks  as  fol- 
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lows 

U('d\r,Q)  =  U(r,Q)-U0)(r,Q), 
where  U^\r,Q)  is  the  incident  field  potential: 


UO)(r,0) 


OO 

^  ^ n  P,n —I  /  2 

n= 1 


(COS  6) 


K(p)^(a.)J 


r  >1, 
r  <  1. 


The  radiation  power  of  horn  antenna  (dipole  radiation  resistor)  was  analysed  using  the  formula 


3^_y _ Kf _ , 

8/fAr2  h  zn(z„2  -0.25)[PZi_1/2(cosy)]2  |  Kz*(sc)  |2  ’ 


(12) 

(13) 

(14) 


c)  d) 

Fig.  2.  Far  field  pattern  for  conical  horn  antenna  at  kl  -  12.5  :  1  -  total  field,  2  -  diffracted  field,  3  -  field  in 

the  case  of  semi  infinite  cone  (a)  y  =  50° ,  kc- 18  ;  b)  y  =  50° ,  kc  -  30 ;  c)  y-  70° ,  kc  =  1 8 ;  d)  y  -  70° , 
kc-  35). 
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where  W  ,  W0  -  (4/3)/rP02k^Z  are  powers,  radiated  by  dipole  in  the  presence  of  conical  horn  antenna  and 
without  it,  respectively. 


3.  The  results  of  the  numerical  analysis 
3.1  Radiated  field  pattern 

Plots  in  Fig.  2  illustrate  converging  forms  of  far  field  distribution  in  the  radiation  zone  of  finite  and 
semi-infinite  cones  in  lighted  angular  sector  with  the  rise  in  the  generatrix  length  of  the  finite  cone.  For 
calculations  we  assumed  that  Pk= l/(4zr)  [A],  Z  =  1  [Ohm],  Forms  of  these  fields,  as  shown  by  studies, 

are  defined  by  the  number  of  modes,  excited  in  semi-infinite  cone.  For  example,  the  field  distribution  of 
semi-infinite  cone,  that  correspond  to  curves  3  of  Fig.  2a,  b,  is  mainly  contributed  by  three  modes  (second 

-  fourth).  The  second  peak  is  formed  because  of  the  third  mode  and  a  final  correlation  between  amplitudes 
of  the  first  and  the  second  peaks  is  established  owing  to  the  fourth  mode.  In  conical  horn  of  finite  length, 
which  far  field  pattern  corresponds  to  curve  1  in  Fig.  2a,  the  third  and  fourth  modes  of  semi-infinite  cone 

are  not  formed  completely  yet.  They  appear 
with  the  rise  in  the  length  of  finite  conical  hom. 
The  second  peak  is  formed  only  at  kcx  >  28  and 
further  rise  in  hom  length  just  weakly  influ¬ 
ences  the  shape  of  far  field  patterns  in  lighted 
sector.  The  rise  in  cones’  aperture  gives  the  in¬ 
crease  of  number  of  modes,  excited  in  semi¬ 
infinite  cone.  So,  the  formation  of  the  same 
field  distribution  by  finite  cone  it  is  necessary 
to  increase  its  length.  This  is  confirmed  by  plots 
on  Fig.  2c,  d.  Therefore,  the  use  of  fields,  scat¬ 
tered  by  semi-infinite  cones  for  analysis  of  far 
field  patterns  for  finite  hom  antennae  in  the  ra¬ 
diation  sector  is  possible,  when  hom  parameters 
9  and  excitation  conditions  provide  the  excitation 

of  modes  that  form  field  distribution  for  corre¬ 
sponding  semi-infinite  cones.  This  property,  in 

Fig.  3.  Influence  of  hom  generatrix  length  on  direc-  its  tum’  allows  re§ulate  field  mode  composition 

.  ..  ,  _  .  of  radiation  from  finite  conical  hom  antennae 

tional  diagrams  at  y,  =50  ,  kl  =  12,5:  I  -  kc,  =17:2  ...  .  ,  .  ..  .  ,  T 

°  ' 1  1  with  using  the  changes  in  its  length.  Let  us  try, 

-  kc\  =15;  3  -  kt\  =13;  4  -  field  distribution  for  with  the  help  of  finite  conical  hom  antenna,  to 

second  mode  of  semi-infinite  cone.  form  the  radiation  distribution  for  the  field  with 

suppressed  higher  modes.  In  this  respect  we  will  study  the  change  of  shape  for  far  field  patterns  from 

conical  hom  antenna  with  the  decrease  of  its  length,  aperture  angle  and  dipole  position  remaining  fixed. 
The  obtained  results  are  shown  in  Fig.  3.  As  can  be  seen  from  plots  in  this  Figure,  the  decrease  in  the  cone 
length  leads  to  approximation  of  shapes  of  the  main  lobes  on  far  field  pattern  of  finite  conical  hom  an¬ 
tenna  and  the  second  mode  of  the  semi-infinite  cone.  Flowever,  they  still  have  substantial  difference  in 
amplitude,  because  of  proximity  of  radiation  source  to  the  aperture. 


3.2.  Near  field  of  conical  horn  antenna 


For  many  technical  problems  the  radiation  characteristics  of  hom  antennas  are  found  basing  on 
Fluygens  method.  The  application  of  this  approach  to  finite  conical  horns  is  grounded  on  the  assumption 
that  the  field  in  hom  apertures  can  be  approximated  by  the  field  of  semi-infinite  cone.  Applicability  of 
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such  approximation  was  studied  in  the  case  of  distribution  of  magnetic  component  of  total  field  on  spheri¬ 
cal  surfaces  with  a  radius  close  to  the  length  of  the  horn.  The  field  calculations  results,  obtained  from  the 
solution  of  wave  diffraction  problems  by  finite  and  semi-infinite  cones  were  analysed.  Their  comparison 
showed  that  upon  approximation  of  source  to  the  vertex  the  total  field  in  aperture  of  finite  horn  antenna 
(i.e.  on  spherical  segment  embracing  cone  basis)  became  closer  to  the  similar  spherical  segment  of  semi- 
infinite  cone  (safe  in  the  region  close  to  the  edge).  This  is  confirmed,  for  example,  by  plots  in  Fig.  4a, 
where  even  for  substantially  wide  horn  we  see  only  small  field  oscillations  in  the  aperture  as  related  to  the 
curve,  corresponding  to  field  distribution  in  semi-infinite  cone.  These  oscillations  are  stipulated  by  wave 
diffraction  on  cone  edge,  and  their  amplitude  increases  with  approximation  to  aperture  boundary.  The  shift 
of  the  source  toward  the  basis  of  horn  leads  to  substantial  difference  in  behavior  of  these  fields,  which  is 
clearly  seen  in  Fig.4b. 


a)  b) 

Fig.  4.  Distribution  of  H  (p  -  field  components  in  vicinity  of  conical  horn  on  a  sphere  of  radius  kr  :  1  - 
whole  field,  2  -  primary  field,  3  -  field  in  semi-infinite  come  (a)  kcx  =  24  ;  kl  -  1,5  ;  yx  =  70° ;  kr  =  24,2 ;  b) 
kcx  =  12  ;  kl  =  10,5 ;  yx  =  20° ;  kr  =  12,2  . 


3.3.  Field  penetration  into  shadow  region 

Using  the  horn  antennae  as  measuring  elements  it  is  important  to  study  the  effects,  stipulated  by 
field  penetration  into  the  shadow  region.  In  order  to  study  them  we  analyzed  the  dependence  of  far  field 
patterns  as  a  function  of  geometrical  parameters  of  the  finite  cone  and  position  of  the  source  under  the 
fixed  values  of  angular  parameter  9  .  Fig.  5a  shows  the  character  of  oscillation  and  decay  of  field  in  the 
direction  9  -  n  I 2  under  the  increase  of  the  horn  length  for  different  cone  aperture  angles,  caused  by 
wave  diffraction  on  the  edge  of  aperture.  It  is  obvious  that  the  screening  effect  should  be  increased  with 
the  decrease  of  angle  near  the  vertex,  which  is  confirmed  by  curves  on  Fig.  5a.  Screening  properties  of 
conical  horn  antenna  in  other  directions  9  are  illustrated  on  Fig.  5b,  where  dependences  of  the  field  in 
certain  directions  on  the  position  of  radiation  source  on  the  axis  are  shown.  As  can  be  seen  from  this  Fig- 
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ure,  the  field  level  and  span  of  oscillations  extinguish  with  transfer  of  observation  point  into  the  shadow 


Fig.  5.  Dependence  of  D(0 )  at  fixed  angle  values  0  =  0t  on  geometric  parameters  of  cone  and  source  posi 
tions;  a)  9-nl 2  ,  (1-  yx  =  70° ;  2—yx  =  59° ;  3-  y,  =  50° );  b)  kc,  =18,  yx= 50°  (1-  di  -  70° ;  2-  9j  =  90° 
3-9/  =110°;  4-  9/  =  150°). 
region. 


Fig.  6.  Dipole  radiation  resistance  in  conical  horn  at  fj=20°;  1  -  kl  =  4,7  ;  2  - 
kl  =  6  ;  3  -  kl  =  7 . 
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Important  characteristic  of  antennae,  which  gives  an  idea  about  their  efficiency,  is  normalized  radia¬ 
tion  resistance.  Its  dependence  on  wave  length  of  horn  at  different  dipole  positions  on  the  axis  is  given  on 
Fig.6.  Curves  on  this  figure  were  calculated  using  formula  (14).  As  can  be  seen  from  this  Figure,  the  level 
of  oscillations  decreases  with  parting  of  radiated  source  from  aperture  into  the  depth  of  the  cone.  The 
dashed  lines  on  Fig.  6  correspond  to  normalized  radiation  resistance  of  radial  electric  dipole  with  kl  -6 
and  kl  -1  on  the  axis  of  semi-infinite  cone  with 
aperture  angle  y,  =  20° .  Therefore,  oscillations  of 
normalized  radiation  resistance  on  the  axis  of  finite 
conical  horn  proceed  in  relation  to  the  value  of  nor¬ 
malized  radiation  resistance  of  dipole  on  the  axis  of 
semi-infinite  cone  with  the  same  aperture  angle  and 
extinguish  with  the  rise  in  the  length  of  the  horn. 

This,  in  particular,  explains  the  possibility  of  ap¬ 
proximation  of  the  field  patterns  of  semi-infinite  and 
finite  horn  antennae  with  the  rise  of  the  length  of  the 
latter. 

Fig.6  shows  that  the  resistance  of  radiation  R  has  its 
distinct  maximum  at  kc-  9,96  and  tends  to  one  if 
kl  >  kc  (this  corresponds  to  the  outlet  of  the  radia¬ 
tion  source  outside  of  horn  antenna).  In  Fig.7  we  see 
the  far  field  patterns  of  finite  conical  horns.  The 
wave  length  of  their  generatrix  takes  a  value  that 
ensures  maximal  radiation  resistance,  and  also  two 
closest  minimums  of  radiation  resistance  of  dipole  in 

cone  at  kl  =  7 ,  y  =  20°  (see  curve  3,  Fig  6).  As  seen 
from  Fig.  7,  the  curve  1  corresponds  to  maximal  normalized  radiation  resistance  and  envelops  the  far  field 
pattern  of  cones,  that  form  minimal  reflection  resistance. 

4.  Conclusions 

1.  The  possibility  of  approximation  of  the  reflection  field  of  the  finite  conical  horn  in  illuminated  area 
by  the  field  of  semi-finite  perfectly  conducting  cone,  is  shown  using  the  example  of  axially  symmetric 
excitation.  The  conditions  on  geometrical  parameters  of  horn  and  the  position  of  the  effective  radiation 
source,  under  which  the  improvement  of  the  approximation  may  be  reached,  are  formulated. 

2.  The  physical  principles  of  the  formation  of  field  in  the  horn  antenna  are  formulated.  It  is  shown 
that  the  best  approximation  is  reached  when  dominant  modes  of  the  semi-infinite  cone  are  formed  in  finite 
conical  horn. 

3.  The  losses  of  the  field  on  radiation  into  the  shadow  area  and  the  energetic  characteristics  of  radia¬ 
tion  of  finite  conical  horn  are  evaluated. 

These  results  will  be  taken  as  a  basis  for  modeling  of  the  behavior  of  the  field  in  the  diagnostic  sys¬ 
tem  under  elaboration. 

5.  References 

1.  A.Z.  Fradin  “Antennae  for  high  frequency  waves”.  Moscow.  Sov.  Radio,  1957,  647p. 

2.  Yu.  P.  Kelohlu,  Yu.  H.  Kozulin,  S.T.  Maclovich  Journal  of  technical  Physics,  1978,  48,  N°5,  p. 

894-904. 


Fig.  7.  Far  field  patterns  of  finite  conical  horn 
antenna  at  kl  =  7  ,  y  =  20° :  1  -  kc-  9.96 ;  2  - 
Ac  =  11.8;  3  -  kc  =  8. 


2/9/2010 


STCU  PROJECT  P-281  -FINAL  REPORT  FF  PAGE  13 

3.  D.B.  Kuryliak,  Z.T.  Nazarchuk  “Analytical  -  Numerical  Methods  in  the  Theory  of  Wave  Diffrac¬ 
tion  on  the  Conical  and  Wedge-Shape  Surfaces.  Kyiv,  Naukova  Dumka,  2006,  275p. 

4.  D.B.  Kuryliak  Izvestiya  Vuzov.  Radioelectronica.  1992,  35,  N°3,  p.  3-9. 


2/9/2010 


STCU 


PROJECT  P-281  -FINAL  REPORT 


PAGE  14 


II.  ESTIMATION  OF  THE  FIELD  MODES  AND  FINDING  THE  INFLUENCE  ON  THE  AN¬ 
TENNA’S  EDGE  AND  DIELECTRIC  INFLUENCE  FOR  THE  FIELD  IN  THE  ANTENNA’S  AP¬ 
ERTURE 


1.  Introduction 

The  wide  use  of  dielectric  and  composite  materials  and  products  requires  creation  of  the  methods 
for  their  non-destructive  testing  (NDT)  and  diagnostics,  particularly  for  determination  of  electrophysical 
and  geometric  parameters.  Classic  methods  of  radio-wave  NDT  allow  performing  defectoscopy  (flaw  de¬ 
tection)  and  determining  one-two  parameters  of  materials  in  the  narrow  range  of  changes  in  these  parame¬ 
ters.  However,  modem  dielectric  materials  and  products  have  a  complex  geometry  and  intrinsic  structure. 
This  complicates  mathematic  models  and  requires  solving  inverse  problems  of  diagnostics,  which,  in  its 
turn,  puts  forward  new  requirements  for  measuring  systems  and  testing  devices.  One  of  the  most  informa¬ 
tive  parameters  in  radio-wave  non-destructive  testing  is  a  complex  reflection  coefficient  (CRC),  since  with 
the  help  of  the  given  frequency  dependence  of  the  CRC  of  electromagnetic  wave,  which  falls  on  the  lay¬ 
ered  dielectric,  the  electrophysical  parameters  and  thicknesses  of  its  layers  can  be  determined. 

During  the  remote  probing  for  the  description  of  the  interaction  between  electromagnetic  field  and 
the  dielectric  medium  the  model  of  the  plane  wave  is  widely  used.  This  model  allows  building  simple  solu¬ 
tions  of  direct  problems  of  probing,  which,  in  its  turn,  strongly  simplifies  the  solution  of  the  inverse  prob¬ 
lem.  However,  the  wave  can  be  considered  as  a  plane  one  only  at  substantial  distance  from  the  source  (far- 
field  region).  In  addition  at  large  distances,  the  signal  reflected  from  the  tested  object  is  weak.  In  order  to 
solve  the  problem  of  measuring  of  CRC  of  layered  dielectric  structure  and  to  determine  their  electrophysi¬ 
cal  parameters  and  thicknesses  of  layers,  the  measuring  scheme  on  the  basis  of  a  doubled  T-bridge  and  its 
mathematical  model  based  on  scattering  matrices  of  independent  modes  have  been  used.  The  work  is 
aimed  at  experimental  determination  of  parameters  for  the  model  of  the  bench  and  study  of  their  influence 
on  measurements  results. 

During  reporting  period  the  following  problems  have  been  solved: 

1.  The  method  and  software  for  automated  calibration  of  the  bench  by  frequency  and  extinction 
have  been  developed.  The  experiments  on  measuring  of  the  length  of  electromagnetic  wave  in  waveguide 
and  calibration  plots  of  the  B-1182  wave-meter  have  been  built,  which  allowed  controlling  the  frequency 
with  the  error  not  exceeding  20  MHz  in  the  53...  78  GHz  range.  The  extinction  of  the  standard  /1.3-38  at¬ 
tenuator  in  the  0...30  dB  range  has  been  measured  and  the  calibration  characteristics  by  extinction  for  all 
the  frequency  range  has  been  built.  The  error  in  extinction  measurements  does  not  exceed  2%. 

2.  The  algorithms  and  programs  for  pre-processing  of  measurements  data  have  been  developed. 
The  spectrum  of  the  digitalized  Ux  and  Uy  signals  of  panoramic  extinction  gauge  has  been  analyzed.  The 
algorithm  for  filtering  and  elimination  of  accidental  and  industrial  interferences  has  been  elaborated.  The 
optimum  step  of  frequency  sampling  has  been  chosen,  the  decimation  (resampling  down)  algorithm  of  fre¬ 
quency  sampling  has  been  developed,  which  allowed,  without  decreasing  the  self-descriptiveness  of  meas¬ 
urements,  shorten  the  volume  of  calculations  on  the  stage  of  determination  of  extinction  and  CRC. 

3.  The  experiments  on  the  CRC  determination  for  standard  loadings  have  been  carried  out.  The 
methodical  and  instrumental  errors  appearing  on  each  stage  of  measurements  have  been  analysed.  On  the 
basis  of  analysis  of  measurements’  results  a  new  mathematic  model  of  interferometer  has  been  proposed. 
This  model  deeper  describes  the  properties  of  the  bridge  scheme,  which  allowed  decreasing  methodical 
errors  of  CRC  measurement.  The  proposed  method  allows  determining  the  module  of  the  CRC  with  an 
error  of  2%  and  a  phase  with  an  accuracy  of  2°  in  5-mm  range  of  electromagnetic  (e.m.)  waves. 

4.  The  calibration  methods  for  the  gauge  system  and  for  measuring  CRC  accounting  the  radiating- 
receiving  aerial  have  been  reviewed.  The  experiments  on  measuring  of  CRCs  for  “horn  antenna  -  plane¬ 
layered  medium”  loading  in  the  62...  65  GHz  range  for  estimation  of  horn  antennas  parameters,  have  been 
performed.  Such  measurements  were  conducted  for  free  space,  metal  plate  screens  and  dielectric  plates  at 
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different  distances  from  antenna  aperture. 

5.  The  algorithm  for  determination  of  scattering  matrix  of  the  radiating-receiving  aerial  on  the  ba¬ 
sis  of  the  measured  CRCs  of  standard  loadings,  situated  in  free  space  is  elaborated.  The  e.m.  wave  in 
plane-layered  medium  is  modelled  as  a  plane  wave  with  normal  incidence  on  the  interface. 

6.  The  full-scale  and  numerical  experiments  have  been  performed,  scheme  models  for  the  aerial  in 
“near-field”  and  “far-field”  regions  have  been  obtained.  The  measured  CRCs  of  standard  loadings  have 
been  compared  with  the  theoretic  ones.  The  dependence  of  modelling  errors  on  the  distance  from  the  sam¬ 
ple  to  the  aperture  of  the  aerial  has  been  analyzed. 


2.  Structure  and  operation  principles  of  the  measuring  bench. 

Creation  of  new  methods  and  means  of  radio-wave  testing  of  multilayer  dielectric  materials  with 
the  use  of  inverse  problems  requires  conducting  precise  experiments  for  verification  of  theoretical  models 
of  electromagnetic  waves  diffraction  on  the  tested  structures  [2,3].  This  needs  the  automated  high  accuracy 
measurement  of  CRC  of  electromagnetic  field  at  multiple  frequencies.  From  this  point  of  view  the  bridge 
methods  of  the  CRC  measurement  have  a  number  of  advantages:  relative  simplicity  of  realisation  of  meas¬ 
uring  scheme,  calibration  simplicity,  wide  frequency  range,  high  measurements’  accuracy,  which  depends 
only  from  the  accuracy  of  calibrated  loading  [5,6]. 

The  laboratory  automated  bench  for  multifrequency  measurements  of  scalar  (voltage  standing 
wave  ratio  (VSWR),  extinction)  and  complex  parameters  (modulus  and  phase  of  reflection  coefficient)  has 
been  created  at  the  PhMI.  This  part  of  the  work  is  aimed  at  refinement  of  the  method  of  measurements  data 
processing,  minimization  of  errors  in  determination  of  CRC  from  the  plane-layered  structure  (PLS). 

The  appearance  of  the  automated  test  bench  for  radio  waves  diagnostics  is  shown  in  Fig.  1. 


The  bench  has  been  built  on  the  basis  of  commercial  panoramic  gauge  of  voltage  standing  wave 
ratio  (VSWR)  and  attenuation  P2-69,  which  main  blocks  are  sweep-frequency  generator  (2)  (see  Fig.  2), 
indicator  unit  (1)  and  SF1F  tract.  Generator  unit  produces  SF1F  oscillations  of  53, 57. ..78, 33  GFlz  frequen¬ 
cies.  The  measurement  of  VSWR  and  extinction  foresees  separate  measurement  of  power  of  incident  wave 
and  a  wave,  which  reflected  or  passed  though  the  structure.  Measurement  of  SF1F  oscillations  power  in  the 
tract  is  performed  with  the  help  of  two  directed  bifurcators  with  detector  sections  (6,7).  The  detected  signal 
of  incident  Uinc  and  reflected  Uref  waves  arrives  at  indicator  unit  R2P-67  (1),  which  forms  the  signal  X, 
proportional  to  the  frequency  of  the  e.m.  wave  and  Y  ~  Uref/Ujnc.  The  indication  and  processing  of  informa¬ 
tion  is  performed  by  a  computer  (4),  joined  though  the  interface  unit  (IU)  (3)  with  X  and  Y  outputs  of  the 
indicator  unit. 

For  the  measurement  of  CRC  the  double  T-bridge  (8)  is  included  into  waveguide  tract  by  the 
“passing”  scheme.  One  of  shoulders  of  the  bridge  is  joined  to  the  unknown  loading  (sample  under  testing) 
(10),  and  the  other  -  to  calibrated  loading  (5).  The  output  signal  of  the  bridge  -  “disbalance  signal”  D  -  is 
a  mixture  of  a  signal,  reflected  from  the  unknown  loading  and  a  reference  wave,  reflected  from  the  cali- 
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brated  loading.  The  calibrated  loading  is  cascade  of  polarization  attenuator  and  phase  changer.  The  pa¬ 
rameters  of  the  model  of  the  bridge  scheme  are  determined  with  the  use  of  standard  loading  -  direct  read¬ 
ing  phase  changer  (movable  short-circuitor  with  a  0.01  mm  scale  step). 


Fig.  2.  Automated  bench  for  multifrequency  probing,  flow  chart: 

1  -  indicator  unit;  2  -  sweep-frequency  generator;  3  -  interface  unit; 

4  -  computer;  5  -  calibrated  loading  with  module  (a)  and  phase  (cp)  con¬ 
trol;  6,1  -  directed  detector  sections;  8  -  double  T-bridge;  9  -  horn-lens 
aerial;  1 0  -  sample;  1 1  -  mechanism  for  moving  the  sample. 

3.  Calibration  of  the  bench  by  frequency  and  extinction 

For  frequency  measurement  of  generated  oscillation  the  waveguide  tract  contains  the  B-1182 
wave-meter.  This  device  has  the  extinction  coefficient  of  about  0. 1  dB  in  all  frequency  range,  and  absorbs 
the  energy  of  electromagnetic  field  in  the  band  <10  MFlz.  On  the  frequency  dependence  of  the  transmis¬ 
sion  coefficient  (TC)  wave-meter  there  is  a  “mark”,  which  position  depends  on  the  wave-meter’s  scale 
(Fig.  3).  The  scale  allows  controlling  frequency  with  the  ~20  MFlz  accuracy.  The  software  of  the  bench 
allows  binding  the  measured  signal  X  to  the  readings  of  the  wave-meter  scale  and  to  frequency. 


dB 


Fig.  3.  Frequency  dependences  of  the  B1 1-82  wave-  Fig.  4.  Scheme  of  wave-meter  calibration:  9  - 
meter,  s  =  3.50;  3.00;  2.50;  2.00  wave-meter;  10  -  short-circuitor;  1 1  -  direct  read¬ 

ing  phase  changer 

In  order  to  build  the  wave-meter  calibration  plot  f  [GFIz]  =  f(s  [mm])  we  have  used  the  bridge 
scheme  (Fig.  4).  One  of  shoulders  of  the  bridge  was  joined  to  short-circuitor,  and  the  other  -  to  direct  read¬ 
ing  phase  changer  (movable  short-circuitor).  By  changing  position  of  movable  short-circuitor  and  balanc- 
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ing  the  scheme  the  wavelength  in  waveguide  was  measured  at  fixed  frequencies.  The  relative  error  of  fre¬ 
quency  determination  in  the  middle  of  frequency  range  is  ~0,07%  (A/.  =  0,005  mm  at  the  wavelength  in 
the  waveguide  of  A.a  =  3,676  mm).  The  calibration  plot  of  the  wave-meter  obtained  by  this  way  is  shown  in 
Fig.  5. 


Fig.  5.  Calibration  plot  of  the  B1 182  wave-meter. 

Measurement  of  the  extinction  Y~UrCf/Ujnc  is  performed  in  seven  dynamic  wavebands  of  the  indi¬ 
cator  unit  (see.  Fig.  2)  with  the  5  dB  step.  For  wideband  signals  the  logarithmic  measurement  regime  is 
foreseen.  The  software  of  the  bench  allows  automatic  pre-processing  of  realizations  Y (X),  conversion  into 
logarithmic  presentation  A  =  10xlg(Y/Y0),  where  Y0(X)  is  a  beforehand  measured  realization,  correspond¬ 
ing  to  the  extinction  level  of  0  dB. 

The  measuring  tract  was  studied  experimentally  with  the  use  of  the  ,11,3-38  polarization  attenuator 
as  an  extinction  standard.  The  measuring  scheme  is  shown  in  Fig.  6.  Fig.  7  illustrates  the  relative  error  of 
the  measured  TC,  as  well  as  the  error  of  standard  213-38  attenuator  (according  to  device  certificate,  dashed 
line).  On  the  basis  of  these  data  one  can  determine  the  relative  error  and  the  linear  range  of  extinction 
measurement.  The  extinction  calibration  plot  A[dB]=A[Y]  has  been  built  (Fig.  8). 


Fig.  6.  Scheme  of  passing  junction  of  at-  Fig.  7.  Relative  error  of  measured  extinction  coefficient, 
tenuator:  12  -  polarization  attenuator 
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The  software  of  the  bench  allows  us  to  process  the  measurements’  results  and  present  them  in  con¬ 
venient  form.  The  measurement  and  calculation  data  are  stored  in  the  PC  memory  and  are  accessible  for 
visual  observation  as  frequency  and  phase  dependences.  Fig.  9a  illustrates  the  measured  dependence  of 
bridge  disbalancing  signal  (in  dB)  as  a  function  of  the  frequency  (MFlz).  Fig.  9b  shows  the  dependence  of 
disbalancing  signal  on  the  phase  of  the  reference  signal  at  the  same  frequency. 


Fig.  9.  Frequency  and  phase  dependences  of  the  disbalancing  signal  of  the  bridge  scheme. 


4.  Model  of  gauge  on  the  basis  of  double  T-bridge 

Traditional  approach  to  the  measurements  of  the  CRC  by  bridge  method  is  in  the  balancing  of 
bridge  scheme  -  comparison  of  the  unknown  CRC  and  calibrated  loading  in  shoulders  of  the  bridge  (see 
Fig.  2).  The  value  of  reflection  coefficient  (module  a  and  phase  cp)  can  be  found  directly  from  readings  of 
the  calibrated  loading. 

For  the  SF1F  range  devices  the  precision  calibrated  loadings  (polarization  attenuators  and  phase 
changers)  are  mechanical  devices,  that  is  why  the  balancing  of  bridge  scheme  at  each  frequency  is  labour¬ 
consuming,  and  under  full  automation  the  use  of  precise  electromechanical  gears  is  needed. 

The  method  for  measuring  of  CRC  with  minimum  (finite)  number  of  manipulations  with  calibrated 
loading,  aimed  at  automation  and  computer  processing  of  measurements  results  is  proposed. 

This  was  salved  by  building  a  mathematic  model  of  the  device,  joining  the  measured  complex  re¬ 
flection  coefficient  V  and  calibrated  loading  Vc(a,cp)  with  disbalance  signal  of  the  bridge  D=Uref/Uinc: 
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D=D(V,Vc(a,cp),S).  (1) 

Here  S  =  S(f)  is  a  complex  scattering  matrix  of  the  double  T-junction,  dependent  on  frequency. 
The  double  T-junction  is  considered  as  a  multiport.  Shoulders  1  and  4  (see  Fig.  10)  are  loaded  on  the  coor¬ 
dinated  units  -  detector  of  incident  Uinc  and  reflected  Uref  wave,  respectively.  Shoulder  3  is  loaded  on  the 
unknown  complex  loading  V(f),  shoulder  2  -  on  calibrated  controllable  loading  Vc(a,(p)  (Fig.  10). 

The  T-joint  is  filled  by  dielectric  (air),  it  is  the  reciprocal  device,  that  is  why  the  matrix  S(f)  is 
symmetric.  Accounting  also  the  symmetry  of  shoulder  of  the  T-joint  one  can  show  that  in  this  scheme  the 
complex  amplitudes  of  signals  in  shoulders  1  and  4  are  bound  by  the  relation: 


+  S2iS42 


V  -  v„ 


{\-S,yV)-{\-SiyVc) 


-s^-v-vc 


(2) 


Fig.  10.  Schematic  model  of  the  bench: 
f  -  sweep-frequency  generator; 

Uref  -  detector  if  the  signal  from  the  re¬ 
flected  wave; 

Vc  -  CRC  of  calibrated  loading; 

K2t  -  scattering  matrix  of  double  T- 
bridge; 

V(f)  -  CRC  of  the  loading  under  meas¬ 
urement; 

S(f)  -  antenna; 

Vs(f)  -  CRC  of  plane-layered  environ¬ 
ment  in  a  free  space 


The  detected  tensions  of  the  incident  and  reflected  waves  are  proportional  to  the  signal  power  in 
the  SHF  tract  U  ~  P  ~  |A|2.  Then  the  normalized  signal  of  bridge  disbalance  can  be  written  as: 


D 


,Vc(a,(p\k(),kl,k2,ki)=  — 


U,-„ 


kr\  +  k\ 


V  -Vn 


(1  -  &3  •  F)  ■  (1  -  &3  -  V)-k2  -  V  -V 


(3) 


where  k0(f),..  k3(f)  are  complex  coefficients,  dependent  on  the  geometry  of  the  2T-bridge  and  the  fre¬ 
quency.  They  describe  properties  of  interferometer  and  can  be  determined  on  the  basis  of  the  full-scale 
experiment. 


5.  Experiment  and  numerical  modelling 


Calibration  of  measurement  scheme  consists  in  determination  of  bridge  scheme  parameters  k0(f),.. 
k3(f)  for  all  frequencies  of  measurement  range.  The  dependence  of  D  signal  on  the  phase  of  calibrated 
loading  cp  for  all  frequencies  of  the  range  and  the  set  of  standard  loadings  (see  Fig.  1 1)  is  conducted. 

Coefficients  of  mathematic  model  of  the  bridge  are  determined  by  the  following  target  function: 


Nv 

F\  (4  v,4 )=  Yjpi  - D(Kt ,  K iai , ^ \  ))’ 


i= 1 


(4) 


where  D,  are  measured  values  of  disbalacing  signal  of  the  set  of  calibrated  loadings  Vc(a,(p)  and  the  known 
standard  loading  Vet;  is  the  number  of  fixed  positions  of  the  phase  of  calibrated  loading.  Optimal  coef- 
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ficients  values  for  the  model  are  determined  by  the  minimum  of  target  function: 

(V-  .,4)*=  inf  4) 

The  direct  search  was  used  as  an  optimization  method. 


(5) 


cpxlO,  deg. 


60  62  64  68  70  f,  GHz. 

Fig.  11.  Measured  phase-frequency  dependence  of  the  disbalancing  signal  D  of  standard  loading. 

Fig.  12  shows  experimental  and  modelled  by  formula  (3)  phase  dependences  of  D  for  five  standard 
loadings  on  the  same  frequency,  for  determination  of  k0(f),..  k3(f)  three  standard  loadings  having  been 
used.  As  can  be  seen  from  the  figure,  this  model  accurately  describes  other  phase  dependences  not  used  in 
the  identification  of  model  parameters. 

Dl/2 


Fig.  12.  Phase  dependences  D:  ■- experiment;  - numerical  modelling  D{pet,Vc( a,(p),&o  )■ 


Value  of  phase  of  standard  loading  depends  on  the  choice  of  the  reference  plane  -  zero  position  of 
phase  changer  of  the  calibrated  loading  (p0-  The  position  of  the  standard  short-circuitor  at  ssh  =  0  mm  was 
chosen  as  a  zero  reference  plane.  Distance  between  two  neighbouring  positions  of  the  standard  short- 
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circuitor  was  AsSh  =  0,5  mm.  Phase  of  the  CRC  of  standard  loading  cpsh  was  determined  through  the  wave¬ 
length  in  the  waveguide  V,  by  the  formula 

cp.s/7  =2tt^,  [rad],  (6) 

ka 

The  CRC  of  the  unknown  loading  is  determined  on  the  basis  of  measured  D;(f)  at  different  values 
of  cpi  for  the  frequence  spectrum;  the  sought  CRC  V*(f)  is  determined  for  each  frequency  by  the  minimum 
of  the  target  function  F2: 

F2(r)=X(a -D{y,Vc{a,yi\k{),...,k2f  ,  V*  =  inf  F2(v).  (7) 

i=1  I^H1 

Fig.  13  shows  frequency  dependences  of  bridge  scheme  parameters  k0(f),..  k3(f). 


J ! ! I_ I_ I_ I 


62.5  63  63.5  64  64.5  65  65.5 


Fig.  13.  Frequency  dependence  of  complex  parameters  of  the  model  of  double  T-bridge  k0..k3(f)  and  phase 
of  zero  plane 

In  order  to  verify  this  method  the  CRCs  of  the  standard  loadings  (short-circuitor)  have  been  meas¬ 
ured.  The  average  deviation  of  the  module  and  phase  along  the  range  do  not  exceed  0.02  relative  unit  (r.u.) 
and  0.04  rad.  The  average  deviation  and  dispersion  of  the  optimal  CRC  V*  from  the  standard  value  by 
module  and  phase  in  the  63... 68  GHz  range  are  presented  in  Table  1 


Table  1. 


Parameter 

Value 

\vet\-\v*\, 

relative  unite  (r.u.) 

0.02 

arg(Fe/)-arg(F*),  rad. 

0.04 

a(|  V*  |)  ,  r.u. 

0.001 

a(arg  V*) ,  rad. 

0.003 
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Let  consider  the  antenna  as  a  passing  linear  device  (four-pole)  with  losses  (Fig.  10).  Let  assume 
that  in  the  free  space  there  is  a  plane  electromagnetic  wave  which  falls  perpendicularly  onto  infinite  (M- 
l)-layered  structure  on  the  semi-infinite  medium  with  flat  interfaces,  which  CRC  depends  on  the  layer 
thickness  and  electrophysical  parameters  of  materials  [3,4] 

Vs=Vs(f,h1,..,hM-i,si,..,SM,tg8i,...,tg8M).  (8) 

Connection  between  the  measured  CRC  V(f)  and  the  CRC  of  the  plane-layered  medium  Vs(f)  can 
be  expressed  as  [5] 


V  =  S1  + 


8*2 

1  —  S^Vg 


or  Vo  = 


Si  -V 


S1S3  -So  -StY 


(9) 


where  Si,.., S3  are  complex  coefficients,  determined  by  the  elements  of  the  scattering  matrix  of  the  antenna 
and  can  be  determined  experimentally. 

In  order  to  determine  the  antenna’s  parameters  as  the  standard  loading  in  the  free  space  Vset(f,  x) 
the  flat  metal  screen  has  been  used,  its  position  was  changed  from  the  zero  plate  x=0  to  20  mm  with  the 
Ax=0.2  mm  step.  Geometric  parameters  of  the  used  antenna:  distance  from  the  aperture  to  the  zero  plane  is 
85  mm,  antenna’s  aperture  was  10x10  mm,  length  of  the  antenna  from  the  mouth  to  aperture  was  50  mm. 

Fig.  14  shows  the  measured  CRC  of  the  “antenna  -  metal  screen”  structure  for  different  positions 
of  the  screen  Vet(f,  x). 


Fig.  14.  Measured  CRC  of  the  antenna-metal  screen  system  Vet(f): _ , _ ,  ....  movement  of  the  screen 

by  Ax=0.4  mm  (different  colours  correspond  to  different  positions  of  the  screen) 


Complex  coefficients  Si(f),..,S3(f)  for  each  frequency  are  determined  with  the  use  of  experimen¬ 
tally  measured  values  Vet(f,  x)  and  theoretic  value  of  the  CRC  of  the  standard  Vsetth(f,  x). 


(Su 


S2,S3) 


=  ^3^, 5^3); 


Nx 

fj(W)=IW-mASi.  vm  (x, )) 

i=0 


where  Xj=iAx.  For  metal  screen  it  is  assumed  that 

| V Set th|  =  1,  arg(Vsetth)  =  I-Itvc/Xq  +  n, 


(10) 


(11) 
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where  x  is  the  position  of  the  screen  compared  to  the  zero  plane,  X0  length  of  the  electromagnetic  wave  in 
vacuum. 

Fig.  15  and  16  show  frequency  dependences  for  coefficients  of  the  antenna  model  Si(f),..,S3(f)  and 
CRC  of  the  standard  loading  recalculated  to  the  free  space  according  to  the  aerial  model  Vs(f)  (9). 


Fig.  15.  Frequency  dependence  of  parameters  of 
the  horn  antenna  model  Si..3(f) 


1.5 


O) 

c 

ra-2 


-4;J — ■ - ■ - ■ - ■ - ■ - ■ - 1 

62.5  63  63.5  64  64.5  65  65.5 

f,  GHz 

Fig.  16.  Measured  CRC  of  the  metal  screen  recal¬ 
culated  for  free  space  Vset(f),  screen  movement  by 
x=0.4  mm 


Fig.  17  shows  the  measured  and  calculated  by  the  model  of  the  plane  wave  [4]  CRC  of  the  fluoro- 
plastic  plate  with  the  following  dimensions:  h=8.0  mm;  s=2.07;  tg§=0.4xl0’ 3.  As  can  be  seen  from  Fig.  16 
and  17,  there  exists  the  correlation  between  theoretical  coefficient  of  reflection  from  the  plane-layered 
structure  and  the  calculated  CRC  accounting  the  influence  of  the  antenna,  the  error  being  about  0.05  r.u. 
The  measurement  error  can  be  stipulated  by  the  error  in  modelling  of  measuring  system,  inadequacy  of  the 
plane  wave  approximation.  The  CRC  in  the  free  space  can  be  affected  by  spurious  reflections  from  the  sur¬ 
rounding  bodies.  As  can  be  seen  from  the  character  of  the  frequency  dependence  Vs(t),  the  reliable  inter¬ 
pretation  of  measurements’  data  according  to  the  plane  wave  model,  one  should  consider  the  maximum 
possible  number  of  frequency  points. 


Fig.  17.  CRC  of  fluoroplastic  plate: 

_  experiment;  ....  the  CRC  calculated  for  the  plane  e.m.  wave 
The  dependence  of  the  model  error  on  the  distance  from  the  aperture  of  the  horn  antenna  to  the 
zero  plane  d  has  been  studied,  the  results  being  presented  in  Fig.  18. 
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Fig.  18.  Dependence  of  modelling  errors  (in  r.u.)  on  the  distance  from  the  antenna’s  aperture  to  the  plane 
of  the  sample  d: 

1  -  c>2(V,  Vth)  of  the  measured  CRC  of  the  “antenna-sample”  system,  Vth=V(S1(f),..,S3(f),Vsth); 

2  -  g2(Vs,  Vsth)  for  the  sample  (Al)  in  free  space  Vs; 

3  -  a2(V,  Vb)  power  of  the  measured  signal  from  the  sample  (Al),  Vb  -  from  the  antenna  in  the  free 

space; 

4  -  a2(V,  Vth)/a2(V,  Vb)  deviation  of  the  measured  CRC  V(f),  normalized  to  the  power  of  the 
measured  signal. 

It  can  be  seen  from  the  figure  that: 

1)  with  the  rise  of  the  distance  the  modelling  error  for  the  measured  CRC  V(f)  (see  plot  1) 
cr(V,  Vth)  decreases. 

2)  modelling  error  of  the  CRC  of  the  structure  in  the  free  space  with  the  rise  in  distance  sharply 
decreases  to  d=20  mm,  and  then  remains  almost  unchanged; 

3)  useful  signal  from  the  structure  (plot  3)  decreases  with  the  distance; 

4)  modelling  error  of  the  measured  CRC,  normalized  on  the  signal  (plot  4)  is  very  close  to  the  plot 
2,  which,  in  the  physical  sense,  is  also  normalised  by  energy  CRC  in  a  free  space  Vs,  equal  here  to  1 . 


7.  Conclusions. 

Automated  bench  for  radio-wave  diagnostics  allows  measurements  in  the  56... 78  GFlz  range,  the 
accuracy  of  frequency  determination  is  ±0.02  GFlz,  extinction  coefficient  is  in  the  range  of  0...30  dB  with 
±1%  accuracy,  complex  reflection  coefficient  with  the  accuracy  of  0.01  relative  unite. 

The  error  of  measuring  of  CRC  with  the  use  of  the  plane  wave  model  for  the  two  simplest  plane¬ 
layered  structures  (flat  metal  screen,  dielectric  plate,  see  Fig.  16,  17,  18)  has  been  determined.  This  error 
was  shown  to  be  of  «  20%  in  amplitude.  In  the  case  of  sufficient  frequency  sampling  the  error  of  model¬ 
ling,  as  can  be  seen  from  Fig.  16,  17,  can  be  considered  as  an  interference  and  for  simple  structure  it 
should  influence  the  diagnostics  results  only  weakly. 

In  order  to  increase  the  accuracy  of  measurement  of  frequency  dependence  of  CRC  the  use  of 
more  realistic  models  for  the  interaction  of  electromagnetic  field  with  layered  environment  is  needed,  in 
particular  they  should  take  into  account  the  change  of  the  power  of  the  reflected  field  and  the  influence  of 
the  finite  dimensions  of  aperture. 

These  results  will  be  taken  as  a  basis  for  modeling  of  the  behavior  of  the  field  in  the  diagnostic  sys¬ 
tem  under  elaboration. 
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III.  ADAPTATION  OF  THE  THEORY  OF  INVERSE  PROBLEM  SOLUTION  TO  ESTABLISH 
THE  PARAMETERS  OF  THE  MULTI-LAYERED  STRUCTURE  TAKING  THE  ANTENNA’S 
GEOMETRY  INTO  ACCOUNT 
1.  Introduction 


The  reliability  enhancement  in  diagnostics  of  plane-layered  dielectric  materials  with  the  help  of 
UHF  waveband  electromagnetic  field  requires  the  enhancement  of  the  experimental  measurement  data  on 
the  basis  of  solution  of  corresponding  inverse  problems.  However,  the  elaboration  of  efficient  mathemati¬ 
cal  methods  for  solution  of  inverse  problems  faces  with  the  principle  difficulties  stipulated  by  their  non¬ 
linearity  and  incorrectness.  One  of  the  ways  to  solve  inverse  problems  of  probing  of  layered  materials  is  to 
determine  parameters  of  layered  structure,  for  which  the  objective  function  for  some  fixed  set  of  frequen¬ 
cies  reaches  the  global  minimum.  However,  despite  of  evident  simplicity  of  such  an  approach,  its  effi¬ 
ciency  substantially  depends  on  the  choice  of  the  model,  initial  approximation,  efficient  minimization  algo¬ 
rithm,  preprocessing  of  measurements  data.  In  addition,  the  objective  function  is  multiextremal,  and  ex¬ 
perimental  data  always  have  errors,  which  requires  modification  of  existing  algorithms  of  development  of 
new  ones,  which  would  provide  reliable  estimation  of  structure  parameters  in  solution  of  particular  practi¬ 
cal  problems. 

2.  Direct  problem  of  probing  of  flat-layered  dielectric  materials 

2.1  Choice  of  physical  model 

The  influence  of  interference  on  measurements  results  lead  to  errors  in  estimation  of  parameters  of 
plane-layered  structure.  Let  consider  here  interferences  as  a  generalized  influence  of  all  possible  factors  on 
determination  of  the  sought-for  parameters.  These  could  be  both  incorrectness  in  the  theoretical  model  and 
measurement  errors.  At  this  stage  of  studies  we  remain  within  the  model  of  the  plane  wave  accounting  that 
the  deviation  of  real  field  of  antenna  from  the  field  of  plane  wave  can  be  accounted  by  introduction  of 
some  level  of  interference  into  the  algorithm  for  solution  of  inverse  problem.  Such  an  assumption  is  ac¬ 
ceptable  at  least  for  simple  structures,  which  has  been  experimentally  demonstrated  in  the  previous  report. 

2.2  Initial  correlations  for  calculation  of  reflection  factor 


Let  us  consider  a  layered  structure  S ,  which  contains  n  - 1  uniform  isotropic  dielectric  layers  with 
flat  boundaries  (Fig.  1),  each  of  which  is  characterized  by  a  thickness  hm  ,  dielectric  permeability  em  and 
the  tangent  of  wave  reflection  angle  8m,  ( m  =  I  ...,n  -  I  ).  Here  we  consider  that  magnetic  permeability  is 
equal  of  that  of  vacuum  p0  .  Let  assume  that  E  -  a  polarized  plane  wave  with  the  e~ja>t  time  dependence 
falls  at  the  angle  0  on  the  upper  limit  z  =  0  of  the  structure  S  .  We  should  find  an  expression  for  reflection 
coefficients  on  the  upper  limit  of  the  structure.  For  this  purpose  let  find  expressions  for  vectors  of  electro¬ 
magnetic  field  in  each  of  layers. 

From  the  geometry  of  the  problem  it  follows  that  the  magnetic  field  has  a  single  component 
// =  H™  (x,z)  in  each  layer,  which  is  normal  to  the  plain  of  the  figure,  and  the  electric  field  has  two 

components,  related  to  the  // by  the  following  correlation: 


K=- 


i 


8H” 


JG)£m  dz  -  J  COE„, 

where  tm  =  sm  (1  +  jdm ) .  The  last  equality  appear  from  Maxwell  equations. 
Than  we  obtain  that 


Em  =  — 


1 


sh; 


dx 
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AH™+k2mH™=0, 

where  km  =  co  e mp0 ,  lmkm  >  0  . 

In  addition  to  the  equality,  the  component  of  magnetic  field  satisfy  the  following  boundary  condi¬ 
tions  on  interfaces  of  uniform  media 


and  conditions  of  radiation  on  infinity  [1],  Here  H‘"c  =/foe/^zcose+xsm0^  is  an  incident  field,  zm  =  ^ h t . 


i= 1 

From  these  we  obtain  the  following  recurrence  relation  for  finding  complex  reflection  coefficient 
for  r0  on  the  surface  of  the  structure  S  : 


rm-\ 


B  +  ea"'  r 

1  +  B  V 


where  Pm  ( 1  pm  )/( 1  ^  Pm  )  ?  Pm  ’  am  ^um^m  ?  ^n—1 

m  =  1, - 1  . 


(1) 


i  > 


2  *  2  2 
^  m  ~  k,n  i 


Re«m  >0, 
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Factors  rm  -  can  be  interpreted  as  reflection  factors  on  the  lower  boundary  of  each  layer.  Only  the 
r0  factor  is  accessible  for  measurement. 

As  we  can  see  from  (1),  the  reflection  factor  r0  is  unambiguously  defined  by  some  combinations 
of  physical  and  geometrical  parameters  em,5 m,hm  .  If  the  changes  of  these  parameters  the  values  of  pm 
and  am  remain  constant,  then  the  characteristics  of  the  field  r0  is  not  changes  as  well.  Thus,  at  certain  value 
the  r0  corresponds  to  a  certain  set  of  real  layered  structures  S .  Let  refer  them  as  equivalent  structures  and 
parameters  pm  and  am  -  effective  parameters. 

Let  us  note  that  the  choice  of  effective  parameters  is  ambiguous.  It  is  evident  that  the  problem  of 
uniqueness  in  determination  of  parameters  of  layered  structure  by  the  r0  reflection  factor  can  be  consid¬ 
ered  only  in  relation  to  effective  parameters.  As  a  rule,  at  some  additional  information  on  a  part  of  physical 
parameters  we  can  reckon  on  unambiguous  determination  of  the  other  part  by  effective  parameters. 

Let  consider  in  details  the  case  of  normal  incidence  of  electromagnetic  wave,  i.e.  9  =  0.  Under 
these  conditions  we  have  um=-jkm  and  pm  =kmsm_]/km_]sm  ,  am  =  2jkmhm  . 

The  values 


8  ft  = 

um  ’ 


■'m- 1 


:  pnPmK 


can  be  considered  as  effective  parameters.  Then  am  =  2j<nvm  yj\  +  jt>m 


Pm 


=  e„ 


1  1 

1  +  fim 


.  These  effective 


parameters  are  generalizations  of  parameters  used  in  [2]. 

There  is  possible  another  way  for  introduction  of  effective  parameters,  which  we  shall  use  during 
further  explanations: 


^ ^ ^  Pm  )  ’  Ym  Pm  j ’  Tm  ,  >  (^) 

"’m— 1 

where  h0  is  some  arbitrary  length.  It  can  be  assumed  equal  to  1.  We  have  aIH  =  co  (  /?.,„  -  y,„  j  and 

pm  =  xm  -  n-1 — 'lulA.  Let  us  note  that  Xm  and  y„,  are  dimensional  parameters:  ] ,  [ym  ]  =  sec.  There- 

P-m  Ym 

fore  for  correct  mathematical  problem  definition  let  introduce  dimensionless  in  the  following  way 

km  —  P)()/-/H  ,  y m  —  m  5  (^) 


where  co0  is  some  fixed  circular  frequency.  Then 


Pm  =  t m  *  Y'”  1  am  =  “  Ym  )  ,  (4) 

Pm  Ym 

where  ©=  K)/co0 .  For  simplicity  sake  further  we  shall  omit 

In  other  words,  each  set  of  equivalent  structures,  which  is  unambiguously  described  by  a  vector 

P  ~  {^0’ 1’  YO’-"’  Yn— lfil’"'! U-l  } 

and  a  certain  value  of  ©  correspond  to  a  concrete  value  of  r0  (©) .  In  a  compact  form  this  can  be  written  as 
an  operator  equation  A (©)[/>]  =  r0  (©) .  For  simplicity  sake  further  we  shall  write  this  operator  equation  as 


Ap  =  r0. 

Operator  A  is  called  the  operator  of  a  direct  problem.  Let  note  that  this  operator  is  nonlinear  and 
continuous.  Its  continuity  is  stipulated  by  fact  that  small  variations  in  structure  parameters  correspond  to 
small  variation  in  parameters  of  electromagnetic  field  and,  correspondingly,  of  reflection  factor. 


3.  Inverse  problem  of  probing  of  layered  structures  and  methods  of  its  solution 


The  problem  of  determination  of  effective  parameters  p  of  layered  structures  on  the  basis  of 
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measured  frequency  dependence  of  reflection  factor  r0  belongs  to  the  class  of  inverse  problems  in  mathe¬ 
matical  physics,  which  in  general  case  are  incorrect.  In  mathematical  aspect  it  means  the  solution  of  opera¬ 
tor  equation  (5).  It  has  been  found  that  when  in  the  equation  (5)  at  a  certain  frequency  range  co  the  precise 
values  of  r0  are  specified  (for  example,  in  the  idealized  numerical  experiment),  the  inverse  problem  has  a 
single  solution  [2], 

Since  the  measured  values  of  r0  contain  not  only  useful  information,  but  as  well  interferences, 
then  is  such  a  case  we  can  talk  only  about  the  approximated  in  a  certain  aspect  solution  of  equation  (5), 
which  is  stable  for  small  changes  of  initial  data  r0  .  For  their  construction  one  can  use  the  additional  infor¬ 
mation  on  the  solution  of  equation  (5).  In  our  case  the  additional  information,  which  has  quantitative  char¬ 
acter,  is  determined  by  the  range  of  changes  of  real  thicknesses  of  structure  layers  and  their  electrophysical 
parameters  and,  correspondingly,  effective  parameters  p  . 

For  description  of  methods  for  solution  of  equation  (5)  in  this  case  let  us  introduce  the  following 
symbols.  Let  put  F  for  a  set  of  possible  values  of  effective  parameters  p  ,  and  R  -  for  a  set  of  possible 

dependences  of  reflection  factor  r0  (co) .  Let  assume  that  F  and  R  are  metric  spaces  (a  certain  metric 
specified  on  these  sets),  and  M  is  a  set  of  possible  solution  of  values  of  effective  parameters  p  of  the 
structure  under  investigation  ( M  eF)  and  N  =  AM  is  an  image  of  the  set  M  after  its  mapping  with  the 
help  of  operator  A  . 

Under  condition  of  r0  ( co)  e  /V  there  exists  an  element  for  which 

Pa  =  A~\  ■  (6) 

The  most  widely  used  method  in  calculation  practice  for  finding  pa  is  a  trial-and-error  method.  It 

means  the  following.  Let  build  the  sequence  {p„},  for  which  p R(Apn,r0)  — » 0  under  condition  of  n  — »  co  . 
The  pn  element  is  considered  as  some  approximation  of  the  solution  p0.  In  [3]  it  has  been  shown  that 
Pn  Pa  >  when  M  is  compact. 

Unfortunately  there  are  no  effective  criteria,  which  would  allow  determine  membership  of  the 
r0(co)  element  in  the  set  N  .  As  a  rule,  it  is  assumed  a  priori.  In  practical  problems  of  diagnostics  r0(co) 
can  be  not  a  member  of  the  set  N  .  In  such  a  case  we  cannot  build  the  approximate  solution  based  on  for¬ 
mula  (6),  since  the  A~lr0  can  be  senseless. 

In  this  situation  one  should  use  generalization  for  the  notion  solution  of  equation  (5)  -  quasisolu¬ 
tion.  According  to  [3]  the  quasisolution  of  the  operator  equation  in  the  set  M  is  an  element  pkeM , 
which  minimizes  the  p R(Apk,r0)  functional: 

Pr(aPicPo )  =inf  P r(aP^o)  ■ 

psM 

If  M  is  a  compact  set,  the  quasisolution  exists  for  arbitrary  r0  ( co)  e  R  and  when,  in  addition, 
r0  (co)  e  N ,  then  the  quasisolution  pk  coincides  with  the  approximate  solution  p{)  =  A  lr() .  At  approxima¬ 
tion  of  the  right  side  of  (5)  to  precise  value  (interference  in  a  metrics  F  tends  to  0)  the  quasisolution  pk 
tends  to  precise  solution  p0 . 

It  is  clear  than  in  our  case  the  set  F  coincides  with  the  3 n  -dimensional  space  C3"  ,  and  R  -  with 
the  space  L2  (<orain,comax) ,  if  r0  (co)  is  set  in  a  range  [<omin,comax]  •  However,  real  measurements  of  reflec¬ 
tion  factors  are  performed  on  discrete  mesh  of  frequencies 

(£>l  =  comin  +  (/-l)Aco,  (1  =  1,..., m), 

where  Aco  is  a  sampling  increment  of  frequency  range.  In  such  a  case  the  set  R  coincides  with  complex 
m  -  dimensional  space  Cm  .  Then  the  problem  of  finding  approximate  quasi-solution  is  reduced  to  minimi¬ 
zation  of  functional 
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m  2 

AP)  =  Z  |  r0(™hP)-re{™l)  |  > 

1=1 

where  re(co are  measured  values  of  reflection  factor  of  electromagnetic  field  on  the  upper  layer  of  the 
structure  under  investigation  for  frequency  co  =  co/ .  As  was  mentioned  above,  proceeding  from  natural  lim¬ 
its  for  parameters  of  layered  dielectric,  one  can  define  the  range  of  changes  in  the  p  parameters  as  a  cer¬ 
tain  limited  region  B0  of  the  space  C3”  ,  which  is  a  compact  set  in  C3"  .  Then  the  inverse  problem  of  prob¬ 
ing  of  layered  materials  is  the  finding  of  the  pk  as  it  is 

f(pk)=  inf  f(p).  (7) 

PgB0 

Since  B0  is  compact,  there  exist  a  single  solution  for  the  problem  (7). 

Numerical  experiments  [7]  indicate  that  the  f{p)  function  has  a  complicated  topologic  structure 
with  a  large  number  of  local  minimums  with  value  of  objective  function  close  to  the  global  minimum,  the 
measure  of  attraction  regions  [11]  of  local  minimum  points  is  relatively  small  in  comparison  with  the 
measure  of  the  sought  set.  Despite  of  intensive  studies  [8-12],  up  to  now  there  are  no  efficient  algorithms 
for  finding  global  minimum  of  linear  functions. 

Concerning  the  practical  applications,  it  is  very  important  to  use  systematic  algorithms  that  guar¬ 
antee  finding  of  ad  lib  close  approximations  of  solution  to  the  finite,  preset  number  of  calculations.  The 
main  drawback  of  methods  of  this  class  is  exponential  [6,15,16]  in  relation  to  problem  characteristics,  es¬ 
timation  of  number  of  calculations,  which  make  their  inefficiency  for  large  scale  problems  or  for  seeking 
high  precision  solutions. 

Therefore  in  practice  one  uses  numerous  stochastic  and  heuristic  methods  that  provide  solution  of 
optimization  problem  only  with  certain  probability  or  do  not  give  any  guarantee  of  it.  As  examples  of  such 
methods  we  can  mention  [5]  cluster  modifications  [9,1 1]  of  random  multistart  [6],  which  combine  methods 
of  local  minimization,  some  strategies  for  generation  of  random  starting  points  and  stochastic  criteria  of 
stopping  [9,1 1]. 

Similar  drawbacks  are  typical  as  well  for  other  methods  used  for  solution  of  inverse  problems,  par¬ 
ticularly  for  genetic  algorithms  [6],  algorithms  of  model  quenching,  neural  nets  [7],  etc. 

In  the  algorithm  based  on  combination  of  the  method  of  nets  [1 1,12]  with  methods  of  local  search, 
has  been  considered.  It  can  be  said  that  this  method  finds  the  point  of  global  minimum  only  in  the  case  of 
the  choice  of  rather  small  meshwidth,  which  requires  large  array  of  calculations  even  for  not  many  vari¬ 
ables.  Here  the  net  is  built  without  accounting  the  global  model  of  behavior  of  the  objective  function  [6], 
therefore  the  convergence  of  this  algorithm  cannot  be  guaranteed.  The  other  drawbacks  of  the  method  is 
the  absence  of  criteria  for  choosing  necessary  precision  for  finishing  of  the  work  of  algorithm  and  estima¬ 
tion  of  error  of  the  obtained  solution  in  the  case  of  incorrect  input  data. 

4.  Numerical  solution  of  inverse  problem  of  electromagnetic  probing  of  layered  dielectrics  by 
the  method  of  branches  and  limits 

4,1  Presentation  of  results  of  mathematical  interpretation 

Characteristic  peculiarity  of  mathematical  theory  for  interpretation  of  data  of  radiowave  probing  of 
layered  structures  is  the  provision  of  uniqueness  of  the  choice  of  the  pk  element,  which  created  the  illu¬ 
sion  of  maximization  of  reliability  of  interpretation  of  results.  During  its  choice  in  addition  of  the  available 
(objective)  a  priori  information,  which  is  already  accounted  during  building  of  the  set  B(} ,  on  the  basis  of 
formal  mathematical  considerations,  an  additional  attribute  is  introduced  into  consideration  (subjectively), 
which  allows  making  preferences  to  some  of  elements  from  the  B0 ,  and  call  the  choice  optimal  on  this  ba¬ 
sis.  In  our  case  such  a  criterion  is  a  minimum  of  misalignment  between  measured  and  theoretical  values  of 
reflection  factor.  The  other  option  of  criterion  is  possible;  it  includes  the  maximization  of  correlation  fiinc- 
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tion  between  precise  and  measured  values  of  reflection  factor. 

Such  a  criterion  for  choosing  optimum  element  provides  convergence  of  the  last  to  the  precise  so¬ 
lution  under  condition  of  tending  of  norm  of  interferences  in  the  measured  values  to  zero.  Of  course,  if  the 
convergence  -  as  a  possibility  to  estimate  parameters  of  layered  structures  with  the  preset  possible  accu¬ 
racy  -  would  be  realized  in  practice,  there  would  be  no  doubts  in  the  efficiency  of  the  form  of  interpreta¬ 
tion  result  presentation  in  the  form  of  a  single  element. 

However  the  assumption  on  tending  of  interferences  to  zero  too  much  extent  appropriate  for  solu¬ 
tion  of  purely  mathematical  issues  of  the  theory  of  inverse  problems.  In  practice  we  always  have  interfer¬ 
ence.  Therefore  if  to  present  the  set  of  approximate  solutions,  which  correspond  to  different  levels  of  inter¬ 
ference  as  a  certain  trajectory,  then  its  behavior  (in  relation  to  the  p0  point)  in  vicinity  of  zero  value  of 
interference  norm,  in  practice  we  are  interested  only  in  the  part  of  this  path  than  correspond  to  the  range  of 
really  possible  values  of  the  interference  norm. 

The  problem  has  as  well  the  following  explanation.  An  arbitrary  set  can  be  completely  represented 
by  one  of  its  elements.  This  allows  stating  that  arbitrary  method  for  interpretation  of  radiowave  probing, 
built  on  the  basis  of  seeking  of  a  single  interpretation  option  from  a  set  of  possible  ones,  does  not  com¬ 
pletely  correspond  to  the  main  aim  of  interpretation  theory  -  the  obtaining  of  the  maximum  volume  of  reli¬ 
able  information,  useful  fro  solution  of  a  practical  problem. 

From  the  practical  aspect  presentation  of  interpretation  results  in  terms  of  norm  of  difference  of 
these  elements 

\Pk~Pa\<e-  (8) 

is  more  efficient  since  interpretation  results  and  precise  solution  of  inverse  problems  are  elements  of  the 
same  space.  In  other  words,  in  addition  to  determination  of  approximate  solution  pk  ,  one  should  find  the 
estimation  of  its  error  (the  e  parameter). 

The  apparatus  for  solving  direct  and  inverse  problems  is  known  to  move  from  analytical  methods 
to  flexible  numerical  methods.  Similarly,  in  the  question  of  efficient  estimation  of  reliability  of  interpreta¬ 
tion  results,  the  analytical  methods  should  make  the  way  to  numerical  ones.  The  reliability  problem  of 
presentation  of  results  should  be  solved  by  algorithm  way  as  a  problem  of  computational  mathematics,  in 
the  first  turn  this  concerns  constrained  extremum  problem  in  determination  of  the  e  parameter  on  the  basis 
of  characteristics  of  interference. 

Therefore  for  finding  the  point  of  the  global  minimum  in  the  problem  (7)  the  modification  of  the 
method  of  branches  and  limits  [8-12]  based  on  the  use  of  interval  analysis  [4,8,9]  has  been  used. 

The  advantages  of  this  method,  in  comparison  with  those  described  above,  is  that  it  belongs  to  the 
class  of  systematic  method,  which  allows  estimating  errors  of  the  obtained  solution  in  the  case  of  inaccu¬ 
rate  input  data  because  of  the  use  of  interval  arithmetic  [4]. 

4.2  Interval  function  of  inclusion  of  objective  function  for  inexact  values  of  reflection  factor 

Before  starting  the  explanation  of  algorithm  for  realization  of  this  idea  let  set  some  definitions, 
which  will  be  used  further  on.  A  set 

X  =  {xeRn  |  -co  <dj  <  X;  <  bj  <  +co  ,  i  =  1 ,...,«} , 

where  Rn  is  n-dimensional  Euclid  space,  a,beRn ,  we  shall  refer  to  as  real  n  -dimensional  interval.  Let 
put  /”  for  the  set  of  n  -dimensional  intervals;  let  consider  X[i]  e  / 1 ,  /  e  {1 to  be  the  ith  component  of 
the  interval  X  e  In  .  For  an  arbitrary  interval  T  =  [xhx2]  e  /  ' ,  -co  <  xl  <  x2  <  co  let  mark  T  =  xl,  T  =  x2.  Let 
call  the  set 

Z  =  X  +  iY  =  {x  +  iy  \xe  X,y  e  Y} ,  where  X,Y  e  In  . 
the  complex  n  -dimensional  interval.  Let  put  C”  for  a  set  of  n  -dimensional  complex  intervals. 

Let  tp(x)is  an  arbitrary  limited  function  defined  on  some  interval  X  el" .  The  interval  inclusion 
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function  [8,10]  of  the  function  cp(x)  on  X  we  call  mapping  ®  :  In  — »  P  ,  for  which  <p(x)  6  0(7)  for  arbi¬ 
trary  range  7cl  and  arbitrary  xeY .  Similarly  to  differentiated  function  we  can  designate  interval  in¬ 
clusion  functions  for  partial  derivatives  of  corresponding  order.  As  we  can  see  the  interval  inclusion  func¬ 
tion  is  not  single. 

The  first  stage  in  building  algorithm  for  solution  of  inverse  problem  of  probing  of  layered  dielec¬ 
trics  and  estimation  of  error  of  the  obtained  solution  is  building  of  inclusion  function  for  f{p) .  Let  con¬ 
sider  the  simplest  option  of  its  construction  with  the  help  of  complex  interval  arithmetic  [4].  The  ways  of 
numerical  realization  of  arithmetic  interval  operations  with  real  and  complex  intervals,  as  well  the  con¬ 
struction  of  interval  inclusion  functions  for  elementary  functions  are  extensively  elucidated  in  literature 
[4,8,10],  therefore  we  shall  not  describe  them  in  details. 

Let  assume  that  we  know  the  absolute  measuring  error  of  r0(co/)  5(co/)  =  Re5(co/)  +  i  Im  5(0)/  ),  val¬ 
ues,  where  Re5(coA)  >  0,Im5(coA.)  >  0  and  some  their  approximate  values  re( co/) .  It  is  clear  that  this  yields 
to  r0((ol)eRe((ol)eIlc  ,  where 

Re(co/)  =  [re( co/)  -  8(a ->i),re{(at)  +  5(co,)]  e  Ilc . 

Let  consider,  instead  of  the  hn,en,5n  e  Rl  values  their  interval  analogues  Hn,En,An  e/1  and  per- 
form  sequentially  all  operation  (1)-  (4)  by  the  rules  of  interval  arithmetic.  Then  we  obtain: 

~  2c0q Mmy]Em\iQ  Re(yi  +  —  2coq Hm yjEm[iQ  Im|^l  +  jAm  j ,  Tm  —  ~~  ? 

1 


P  =  T 

M  m  m 


K  =  HJK,  -  ) ,  B  =  (1  -  Pm  )/(l  +  Pm);  Rm_  J  =  A 

JK-K  1  +  B  meA-Rm 


R  -\-e  mR 


)  Rn- 1  _  B„_i  ■ 


Then  interval  inclusion  function  for  objective  function  can  be  written  as 


f(p)  =  Z  |r0(®/;^)-^(®/)  |  • 

i=i 

It  allows  obtaining  the  lower  and  upper  estimations  of  values  of  objective  function  (2)  on  the  arbi¬ 
trary  interval  P  c  R0  :  / ( p)  e  ^F( P), F(P)J  for  arbitrary  p  eP and  error  6(mk ) . 


4.3  Algorithm  for  solution  of  inverse  problem  of  probing  of  layered  structures 


Let  introduce  some  auxiliary  symbols.  Let  consider  sets  of  elements  in  the  form  e  =  (  B, a,p  ), 
where  Bel3'1  ,a,fieR  are  lower  and  upper  limits  of  estimation  of  the  /( p)  function  at  the  B  interval, 
respectively.  Let  put  L  for  a  set  of  such  elements,  for  which  in  intervals  B  there  is  not  yet  a  decision  on 
possible  location  of  problem  solution  in  them,  and  G  -  for  a  set  of  such  elements,  that  ||i?||  <  A0 ,  where 

0  <  A0  <  1  is  an  arbitrary  constant,  ||i?||  =  maxR[/]  -  B[i]  ,BeIn.  Let  mark  the  upper  estimation  of  the  possi- 

i 

ble  value  of  the  objective  function  on  B0  as  z  .  Algorithm  contains  the  following  steps: 

0.  Set  0  <  A0  <  1,  0  <  q  <  1 .  Assume  G  =  0 .  Calculate  F(B0) ,  assume 
y  =  F(B0),z  =  F(B0),L  =  {(B0,y,z)}. 

1.  If  L  =  0  ,  then  go  to  the  step  7,  otherwise  find  an  arbitrary  element  e  =  (  /i,a,[l  )  within  L  with  the 
minimum  value  of  a .  Assume  L  =  L\{e} . 

2.  If  a  >  z  ,  then  exclude  element  e  from  consideration  and  go  to  step  1. 

3.  If  \\B\\  <  A0 ,  then  assume  G  =  Gu{e}  and  go  to  step  1  1 . 

4.  If  (F(mid  B)-  z)  I  z>  q  ,  then  go  to  step  5,  otherwise  find  some  point  of  local  minimum  x*  eB  of  the 
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function  /  ( p } .  Calculate  /*  =  F(x*)  .  Assume  z  =  min{z,f*} . 

5.  Split  the  interval  B  into  two  -  Bx  and  /i2  by  the  largest  coordinate: 

B  =  By  u7?2,  Bx  nB2  =  0  . 

6.  Calculate  F(By),F(B2) ;  assume  a,  =  T(£() ,  [1,  =  F(Bi) , i  e  {1,2} ;  z  =  min  {z,(31,p2  }.  Assume 
I  =  Iu{(5;,a;,p,)} ,  i  e  {1,2} .  Go  to  step  1. 

7.  Build  such  an  interval  P*  e  73"  ,  that  P*[i]  =  [mine[i],maxe[i]],  i  =  l,...,3n  . 

eeG  eeG 

8.  The  end  of  algorithm. 

The  q  parameters,  used  in  step  4,  serve  for  controlling  of  the  number  of  application  of  the  local 
minimization  method.  The  method  of  local  minimization  is  applied  in  an  interval  B  only  when  the  upper 
limit  of  inclusion  function  in  the  central  point  of  this  interval  is  smaller  or  comparable  with  z. 

Proof  of  convergence  of  the  proposed  algorithm  is  similar  to  the  proof  of  convergence  of  the  inter¬ 
val  method  of  branches  and  limits  [8,10,12]. 

For  analysis  of  the  algorithm  we  have  performed  the  numerical  experiment  for  determination  of 
thickness  of  layers  in  dielectric  structures  with  the  following  parameters:  a)  h}  =  0.02  ,  et  =  3.9  ,  5,  =  0.01  ; 
b)  hx  =  0.01,/z2  =0.02,  £j  =3.9,e2  =2.5,  Sj  =52  =0.01.  The  values  of  reflection  factor,  calculated  by  for¬ 
mula  re(co/)  =  r0(co/)  (1  +  55/100),  where  5  is  evenly  distributed  on  [-l,+l]  arbitrary  value,  5  is  relative 
error  in  percent,  in  the  range  of  relative  frequencies  from  comin  =0.5  to  comax  =4.1 ,  which  at  co0  =  420.0 
corresponds  to  the  frequency  range  from  10  to  82  GFlz,  was  assumed  to  be  the  measured  ones. 

During  numerical  experiment  we  studied  the  dependence  of  the  error  of  the  obtained  solution  by 
the  argument  Ax  on  the  number  of  probing  frequencies  N  and  relative  error  8  of  input  data.  We  accepted 
the  following  values  of  parameters:  q  =  10  6  and  A0  =  10~5  .  The  table  provides  the  averaged  values  of  er¬ 
rors  (after  50  trials). 


Structure  a) 

Structure  b) 

n\ 

10 

20 

30 

40 

10 

20 

30 

40 

5,%\ 

0 

9.8E-6 

9.8E-6 

9.8E-6 

9.8E-6 

1.5E-5 

2.2E-5 

2.2E-5 

2.2E-5 

2.5 

3.0E-5 

4.7E-5 

2.9E-5 

3.1E-5 

1.9E-4 

1.8E-4 

1.7E-4 

1.7E-4 

5 

5.9E-5 

7.9E-5 

5.9E-5 

5.9E-5 

3.5E-4 

3.3E-4 

3.2E-4 

3.2E-4 

10 

1.1E-4 

1.5E-4 

1.2E-4 

1.1E-4 

1.8E-2 

6.4E-4 

9.1E-3 

6.5E-4 

15 

1.8E-4 

2.3  E-4 

1.8E-4 

1.8E-4 

2.9E-2 

2.55E-2 

2.32E-2 

2.02E-2 

20 

8.5E-3 

3.1E-4 

2.4E-4 

2.4E-4 

3.4E-2 

3.1E-2 

3.1E-2 

3.1E-2 

The  presented  results  confirm  the  conclusion  on  the  efficiency  of  application  of  interval  method  of 
branches  and  limits  for  practical  solution  of  problems  in  electromagnetic  diagnostics  of  layered  dielectrics 
for  not  many  unknown  parameters  and  insubstantial  errors  in  input  data. 

Preliminary  analysis  of  results  of  numerical  modeling  shows  that  the  field  of  application  of  this 
method  can  be  substantially  widened  (structures  with  larger  number  of  layers,  decrease  in  requirements  to 
the  accuracy  of  input  data)  by  the  following  way. 

First,  this  is  the  application  of  e  inclusion  functions  built  with  the  help  of  centered  forms  [4,12]; 
second,  for  quick  finding  the  value  of  objective  function,  close  to  the  global  minimum,  one  can  use  the 
tunnel  algorithms  [4,8]  based  on  tunnel  functions  (filling  functions). 

Then,  elements  of  the  set  G  obtained  as  a  result  of  the  work  of  algorithm,  are  located  in  vicinity  of 
the  point  of  the  global  minimum,  forming  the  so-called  cluster  [8].  Accounting  large  number  of  such  ele¬ 
ments,  particularly  in  the  case  of  inaccurate  measurement  data,  it  would  be  promising  to  use  the  modified 
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method  of  branches  and  limits  by  the  strategy  on  isolation  of  points  of  local  minimum  [8]. 

5.  Conclusions 


The  theoretic  study  of  direct  and  inverse  problems  of  electromagnetic  probing  of  multilayered 
structures  has  been  performed. 

It  has  been  found  that  on  the  basis  of  frequency  dependence  of  reflection  factor  it  is  possible  to  de¬ 
termine  only  effective  parameters,  than  describe  some  set  of  real  layered  structures  (  equivalent  structures). 

The  problem  of  determination  of  effective  parameters  of  layered  structures  is  reduced  to  the  solu¬ 
tion  of  inverse  problem,  which  requires  the  minimization  of  misalignment  functional  (objective  function) 
measured  and  theoretically  calculated  reflection  factors  on  the  surface  of  structure. 

It  has  been  shown  that  in  the  presence  of  interferences  in  addition  to  determination  of  approximate 
solution  one  should  find  the  estimation  of  its  error.  For  solution  of  this  problem  a  numerical  method  based 
on  a  new  algorithm  for  global  minimization  of  objective  function,  which  is  the  combination  of  interval 
method  of  branches  and  limits  with  methods  of  local  search. 
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IV.  MODIFICATION  OF  THE  SOFTWARE  FOR  THE  SOLUTION  OF  THE  INVERSE  PROB¬ 
LEM 


Branches  and  limits  method  [8-12]  and  its  application  for  study  of  multiextremal  mathematical 
models  have  been  widely  spread  during  recent  decades.  It  is  one  of  a  few  systematic  methods,  applicable 
for  solution  of  a  wide  spectrum  of  continuous  and  discrete  optimization  problems  under  relatively  simple 
assumptions  concerning  an  objective  function  and  a  search-set. 

Principle  idea  of  the  method  is  a  step-by-step  search-set  partition  of  the  initial  problem  into  several 
subsets  with  subsequent  exclusion  of  subsets,  which  do  not  contain  the  global  minimum  and  continuation 
of  search  in  subsets  considered  as  “perspective”  in  scopes  of  possible  finding  in  them  the  solution  under 
certain  criterion. 

The  substantial  drawback  of  the  method  of  branches  and  limits  is  an  exponentional  growth  of  the 
amount  of  calculations  in  relation  of  problem  characteristics.  Let  note,  that  such  a  pessimistic  estimation  is 
typical  of  all  systematic  methods.  However,  if  to  apply  interval  analysis  and  means  to  accelerate  the  con¬ 
vergence  for  the  method  of  branches  and  limits,  numerical  experiments  confirm  the  possibility  of  it  effi¬ 
cient  realization  [4,8,9]  for  optimization  problems. 

Let  enter  for  convenience  sake  some  designations.  The  set 
X={xeRn  |  -co<  a, ■<xi<bi  <  +co  ,  i  =  l,...,n} ,  where  a,beR"  are  some  vectors,  let  refer  to  as  n  - 
dimensional  interval.  The  set  of  all  n  -dimension  intervals  by  In  ,  and  consider  X[i\  e  /  ,  i  e  {1 ,...,«}  as  the 
ith  component  of  the  interval  X.  For  arbitrary  single-dimension  interval  T  =  [j|,x2]  6  /  ,  -oo  <  x1  <  x2  <  <x> 
identify  T_  =  x1 ,  T  =  x2  1 . 

Let  (p(x)  is  an  arbitrary  bounded  function,  defined  on  a  certain  interval  X  e  Then  the  interval 
inclusion  function  of  the  function  ip(x)  on  X  we  consider  such  an  image  0:7"  — >/'  that  (p(x)eO(T) 
for  arbitrary  interval  7cl  and  arbitrary  xeY  .  Similarly,  for  differentiated  function  we  can  denote  in¬ 
terval  inclusion  functions  for  partial  derivatives  of  corresponding  order. 

Interval  inclusion  functions  are  widely  used  in  modifications  of  the  method  of  branches  and  limits  for  solu¬ 
tion  of  optimization  problems,  since  by  definition  they  allow  easy  determining  upper  and  lower  limits  for 
objective  functions  and  its  derivatives. 

There  is  substantial  obstacle  on  the  way  of  efficient  application  of  the  method  of  branches  and  lim¬ 
its  for  solution  of  applied  problems,  which  is  a  cluster  problem  [8].  Its  essence  is  that  during  the  work  of 
the  method  there  appears  the  necessity  to  analyze  a  large  number  of  subset  with  smaller,  compared  to  the 
initial  set,  diameter  (their  aggregation  is  usually  referred  to  as  cluster),  which  are  situated  in  vicinity  of  the 
global  minimum  point.  Here  the  main  part  of  the  calculations  work-time,  performed  by  the  method  of 
branches  and  limits,  is  lost  for  exclusion  of  these  subsets. 

In  the  case  of  successful  solution  of  the  cluster  problem  we  can  hope  to  enhance  the  efficiency  of 
the  method  concerning  most  of  the  applied  problems. 

It  is  clear,  that  there  is  now  necessity  of  finding  the  point  of  the  global  minimum  with  high  accu¬ 
racy.  This  is  stipulated  by  substantial  errors  in  measurement  data. 

Then,  in  order  to  overcome  the  cluster  problem,  we  shall  use  the  algorithm  in  method  of  branches  and  lim¬ 
its,  that  allows  finding  points  of  local  minimum  and  build  some  circles  around  them  which  contain  local 
minimum.  The  method,  in  its  turn,  should  be  modified  in  the  way  that  the  partition  of  the  search  set  would 
occur  mostly  outside  the  built  circles.  This  idea  is  geometrically  illustrated  in  Fig.  1 

Let  f(x)eC2(B0)  ,B0  e/”  and  all  local  minimum  points  of  the  f(x)  function  on  B{]  arc  internal. 
We  should  find  such  a  point  xgi  e  B0  that 

f(xg/)  =  mmf(x). 
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Fig.  1 .  Line  of  level  for  multiextremal  function  and  strategy  for  solution  of  the  cluster  problem. 


Realization  of  the  proposed  above  idea  of  finding  the  point  xgl  needs  development  of  constructive 
and  efficient  method  for  solving  such  problems  for  arbitrary  interval  B  e  B{) : 

1)  find  a  point  x*  e  B  of  the  local  minimum  of  the  f(x)  function  on  B  ; 

2)  build  a  certain  interval  5'cS,  that  x*  eB  and  x*  would  be  a  local  minimum  of  / (x)  on  B' ; 

3)  efficiently  “cut  out”  the  hyperinterval  B'  from  B  ,  i.e.  build  a  set  of  subintervals  Bl,B2,....,Bm  , 
m  >  1  for  which  Bt  c5,  B  =  B'[jBi,  i  =  1 

For  solution  of  the  Problem  1)  one  can  use  any  globally  convergent  algorithm  of  conditional  local 
minimization.  For  efficient  accounting  simple  restriction  x  6  B  it  worth  to  use  quasi-Neuton  modification 
of  the  linearization  method  [1-3]. 

Let  consider  the  Problem  2).  LetXe/"  is  arbitrary  hyperinterval  (interval  vector)  and 
wX  =  max  [Ti-TA.  In  order  to  solve  the  problem  it  is  enough  to  build  the  such  a  hyperinterval  B'  a  B  , 

that  x*  =  midB'  and  wB'  =  A  ,  where  midX  is  the  central  point  of  hyperinterval  X  ,  and  A  is  the  preset  pre¬ 
cision  of  the  solution  by  argument. 

The  algorithm  A1  described  below  allows  to  find  the  solution  of  the  problem  2). 


0.  Define  intervals  B  ,B' .  Assume  5  =  0,  i  =  1 . 

1.  Suppose  T\n  =  \B\i\B'\iM  ,  T[j]  =  B[j]  for  j  =  \,...,n,j . 

Calculate  F(T)  and  suppose  a  =  F_(T),  [3  =  F(T). 

2.  Suppose  B  =  B\  {T} ,  5  =  5U{  (T,a,|3)  } . 

3.  Suppose  T[i]  =  [5'[i];5[i]] ,  T[  j]  =  B[  j]  for  j  =  \,...,n,j  *i .  Calculate  F(T)  and  assume 
a  =  F(T) ,  (3  =  F(T). 

4.  Suppose  B  =  B\{T},  5  =  5U{(T,a,p)}. 

5.  If  i  <  n ,  then  assume  i  =  i  + 1  and  go  to  step  1 . 

6.  End  of  algorithm. 


From  the  building  of  algorithm  it  is  seen  that  as  a  result  of  its  work  the  set  5  will  contain  not  more 
than  2 n  elements,  which  allows  to  assume  it  efficient  for  large  n  . 


1.  Algorithm  for  solution  of  reverse  problem  of  probing  layered  structures 


The  first  stage  in  building  the  algorithm  for  solution  of  reverse  problem  of  probing  layered  dielec¬ 
trics  and  estimation  of  an  error  of  the  found  solution  is  the  construction  of  the  inclusion  function  for  the 
object  function.  Let  consider  the  simplest  case  of  its  building  with  the  help  of  complex  interval  arithmetic 
[4],  The  methods  for  numerical  realization  of  arithmetic  interval  operation  with  real  and  complex  intervals, 
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as  well  as  the  building  of  interval  function  of  inclusion  for  elementary  functions  are  widely  presented  in 
literature  [4],  so  we  shall  not  describe  them  here. 

Let  assume  that  the  absolute  error  of  measurements  of  values  of  r0( co/)  8(co/)  =  ReS(co/)  +  /Im8(co/), 
where  Re5(coA.)  > 0,Im8(coA)  >  0  and  some  their  approximates  /;,(«/)  are  known.  It  is  obvious  that  then 

>i)(o)i  )  e  Re(o)/)  e  f'c  ,  where 

Re(co/)  =  [re(co/)  -  5(co/),re(co/)  +  5(co/)]  e  Ixc . 

Instead  the  hn, e„,5„  eR1  values  let  consider  their  interval  analogues  Hn,En,An  e  /'  and  sequen¬ 
tially  perform  all  operation  in  (1)-  (4)  (see  the  T03  report)  according  to  rules  of  interval  arithmetic.  Then 
we  obtain: 


■4 m  2©q Hm-\J EmHo  Re  |^1  +  jAm  j ,  Mm  2o)^l I nt  I  m  (  y  I  +  jAm  j 


T  =  .H» 


H 


m- 1 


pm  =  T JA'"-1  Am  =  <b(jAm  -  Mm ) ,  Bm  =  (1  - Pm )/( 1  +  Pm);  Rm_x  = 

Then  interval  inclusion  function  for  the  object  function  will  look  like  [5-7] 

m  2 
F(P)  =  Z  |  *o(©,;P) -**(©,)  |  ■ 

/= l 


B  +  eA"'  R 
1  +  B  e A"'  R 


Rfi— 1  ®n-l  • 


It  allows  obtaining  lower  and  upper  estimations  of  the  object  function  (2)  on  arbitrary  interval  P<^B0: 
l(T’),T’(P)J  for  arbitrary  p  e  Rand  error  5(coA) . 

Let  introduce  some  auxiliary  notations.  We  shall  consider  sets  in  the  form  E  =  (B, a,p  )  ,  where 
B  e  Iin ,  a,p  eR  is  the  lower  and  upper  estimations  of  the  value  of  the  function  /  ( p)  on  the  interval  B  , 
respectively.  Let  put  L  for  a  set  of  such  elements,  for  which  on  intervals  B  there  is  still  no  decision  on 
possible  finding  solution  in  them,  and  G  for  a  set  of  intervals,  which  width  is  smaller  than  some  small 
value  8  >  0 

Together  with  these  sets  we  consider  the  set  C  made  of  elements  in  the  E  =  (  ZLaJLx*  )  form,  where 


eR"  is  a  certain  point  of  the  local  minimum  of  empirical  point  function  of  objective  function  f(x) . 


B  e  I" ,  a,p  g  R  have  the  same  sense. 

Let  denote  the  upper  estimation  of  the  minimum  value  of  target  function  on  B0  as  z  .  The  algo¬ 
rithm  has  the  following  steps. 

Description  of  algorithm  parameters. 

A  >  0  is  a  fixed  solution  precision  by  argument. 

0  <  y  <  1  is  a  fragmentation  constant. 

0  <  8  <  1  is  a  limit  of  hyperinterval.  During  the  performance  of  the  algorithm  the  hyperinterval  B  shifts  to 
the  set  G,  if  wB<  A 

0  <  s0  <  1;  Sq  «  8  —  utter  limit  of  partitioning  of  hyperinterval.  The  partitioning  of  hyperintervals  is  forbid¬ 
den  if  wB  <  s0  . 

Let  consider  that  the  problem  is  solved,  if  at  certain  step  their  is  fulfilled  condition 

wP  <  A, 


Where  P  =  (Pl,..,Pn)eIn ,  Pt  = 


min  X_h  max  X  / 

XglUcUg  XeiUcUG 


,  i  = 


Condition  for  the  exclusion  of  sets  from  consideration  takes  the  form:  F(X)  >  0  where  F(  X)  is  any  in¬ 
terval  function  of  inclusion  of  object  function. 

The  algorithm  is  the  anti-cluster  modification  of  interval  method  of  branches  and  limits  of  the  0-order  in- 
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eludes  the  following  steps. 

0.  [initialization].  To  chose  values  of  A  >  0,  e  >  0,  e0  >  0,  0  <  q  <  1,  0  <  y  <  1  parameters;  calculate 
F(B0) ;  assume  y  =  F(B0) ,  z  =  F(B0 ) ,  L  =  {  ( B0,y,z )} ,  G  =  0,  C  =  0. 

1.  [test  to  achieve  accuracy].  If  wP  <  A  ,  then  go  to  step  10. 

2.  [refinement  of  clusters].  Iff  =  0  Ta  C  ^0  then  go  to  step  9. 

3.  [refinement  by  accuracy].  If  L  =  0, C  =  0, G  + 0  and  e >  e0 ,  then  assume  e  =  s  •  y ,L  =  G,G  =  0  and 
go  to  step  1 .  Otherwise,  go  to  step  1 1  . 

4.  [search  of  the  perspective  element].  Find  in  L  an  arbitrary  element  E  =  (B,  a, [3  ),  for  which 
wE  =  max  wX  (element  of  the  maximal  width).  Assume  L=L  \  \E  \ . 

Ie  L 

5.  [tolerance  test]  If  a  >  0  ,  then  exclude  the  element  E  from  further  consideration  and  go  to  step  1 . 

6.  [optimality  test]  If  wB  <  s ,  then  assume  G  =  G  U  {E}  and  go  to  step  1 . 

7.  [modified  partitioning] 

7.1.  [condition  of  modified  partitioning].  If  (F (B)  -  z)  /  z  >  q ,  then  go  to  step  8. 

7.2.  Using  the  method  of  local  minimization,  find  a  point  of  local  minimum  x*  eB  of  the  function 
f(x) .  Calculate  f*  =  /(x*) .  Assume  z  =  min{z,f*} . 

7.3.  Build  an  interval  5'cfi,  x*  =  mid  B',  wB'  =  A. 

7.4.  If  FJB')  >  0  ,  then  exclude  B'  . 

7.5.  Using  algorithm  Al,  build  a  partition  S  of  the  intervals.  Assume  L  =  L{JS  and 
C  =  C{J{B'} 

7.6.  Go  to  step  1. 

8.  [standard  partition].  To  find  index  k  for  which  wBk  =  max  wBj .  Divide  the  hyperinterval  B  into 

two  hyperintervals  by  kth  coordinate:  B  =  Bx  \JB2  .  Calculate  F( Bl ) ,  F( B2 ) ;  assume  a,-  =  F( Bi ) , 
P,-  =F(Bj ) ;  z  =  min  {  z,pl5p2  }  ,L=L[j  {(S^a^p,)} .  Go  to  step  1. 

9.  [treatment  of  the  set  C],  If  C  =  0  ,  then  go  to  step  1 . 

9.1.  Choose  an  arbitrary  element  E  =  (B,  a,  p,  x*  )  e  C  . 

9.2. If  a  >  0  ,  then  go  to  step  9.7. 

9. 3. If  wB  <  e ,  then  assume  G  =  G  U  {E}  and  go  to  step  9.7. 

9.4.  Build  the  interval  5'cfi,  x*  =  mid  B',  wB'  =  wB  *  y. 

9. 5.  Assume  B  =  B'  ?  a  =  F_{B) ,  p  =  F(B),  z  =  min  {z,p1,p2  }  and  go  to  step  9.1. 

9. 6.  Using  the  algorithm  A2,  build  a  partition  S  of  the  interval  B  .  Assume  L  =  L{JS  . 

9.7.  Assume  C  =  C\  {E}  and  go  to  step  9. 

10.  [precision  is  achieved].  Go  to  12. 

1 1 .  [precision  is  not  achieved]. 

12.  [end].  Print-out  the  hyper-rectangle  P.  End  of  algorithm. 

Parameter  q  and  checking  at  the  step  7. 1  serve  for  controlling  the  number  of  the  applied  modified  par¬ 
tition.  Modified  partition  in  the  algorithm  is  used  for  certain  interval  B  only  when  the  upper  limit  of  the 
inclusion  function  on  this  interval  is  smaller  or  commensurable  with  z. 

2.  Program  realization  of  the  modification  of  interval  method  of  branches  and  limits  for  solu¬ 
tion  of  inverse  problem  of  probing  of  layered  structures 

The  described  above  algorithm  is  realized  with  the  help  of  methods  of  object-oriented  program¬ 
ming  in  C++  as  a  batch  of  application  programs.  Structurally  the  batch  contains  two  parts:  auxiliary  pro¬ 
gram  DirTask  and  main  program  DielScan. 
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The  DirTask  program  serves  for  calculation  and  recording  in  the  text  file  values  of  reflection  fac¬ 
tors  and  results  of  modeling  of  errors  of  their  measurements  for  the  preset  structure  and  chosen  frequency 
range  (imitation  of  reflection  factor  measurements  data).  These  data  will  be  input  data  for  the  DielScan 
program.  Graphical  interface  of  the  user  of  the  DirTask  program  is  shown  in  Fig.  2.  It  contains  fields  for 
input  of  parameters  of  dielectric  structure  for  calculation  of  complex  reflection  factor  and  other  auxiliary 
information. 


General  parameters  Discrete  uniform  frequency  net' 

_  'Frequency  range,  GHz 

Problem  _ _ 

mb  |  0.5  ma*  |  4.2 


Number  of  range  divisions  10 


Problem  name 

Number  of  internal 
layers  of  structure 


Parameters  of  method  of  branches  and  limits 
Solution  accuracy  by  argument]  ie-005 
Max  number  of  division 
elements  (Unchecked)  •  u  "" 

Max  number  of  division  [^000000 
elements  (Good) 

Division  limit  | 19-011 


Parameters  of  dielectric  structure 


Errors  of  reflection  factor 


Path  of  file  formation 


Form  files 


>  I  D:\Tes  tProbtems\ 

Value  of 

errors 

5 

T 

Type  of 
errors 

pel 

Fig.2.  General  outlook  of  the  interface  window  of  the  DirTask  program. 


Parameters,  which  values  are  needed  for  the  work  of  this  program  are  divided  into  several  groups. 
General  parameters.  Field  «Problem  name»  contains  the  name  of  text  file,  in  which  the  output  information 
will  be  saved.  The  field  «Number  of  internal  layers  in  the  structure))  gives  the  number  of  layers. 

Discrete  uniform  frequency  net.  The  values  of  lower  and  upper  limits  of  the  chosen  frequency  range  are 
put  into  the  fields  «min»  and  «max)>,  respectively.  The  field  «The  number  of  range  divisions))  contains  the 
number  of  intervals,  into  which  the  frequency  range  should  be  divided. 

Parameters  of  the  method  of  branches  and  limits.  This  group  of  fields  is  provided  by  necessary  in¬ 
formation  for  the  methods  of  branches  and  limits.  The  value  of  A  is  given  in  the  field  ((Solution  accuracy 
by  argument)),  and  the  value  of  the  utter  division  limit  e0  is  put  into  the  field  ((Division  limit  ».  Numbers 
of  maximum  quantities  of  elements  in  L  and  G  sets  are  contained  by  fields  ((Max  quantity  of  division 

elements  (Unchecked)))  and  «Max  quantity  of  division  elements  (Good)»,  respectively. 

Parameters  of  dielectric  structure.  Electrophysical  (dielectric  constant  8  ,  dielectric  dissipation 
tgd )  and  geometrical  (layer  thickness  h )  parameters  for  layers  of  structures  and  errors  of  their  determina¬ 
tion  are  provided  through  fields  of  the  table  in  the  central  part  of  the  window.  In  the  same  table  parameters 
of  structure,  assumed  to  unknown  should  be  marked  by  «?». 

In  the  lower  part  of  the  window  there  should  be  filled  a  field  ((Path  of  file  formation))  with  the 
subdirectory,  in  which  the  files  with  results  of  calculations  by  the  DirTask  program  should  be  stored.  Type 
of  error  in  reflectivity  factor  is  chosen  from  the  floating  list:  «Rel»  or  «Abs»,  and  its  value  is  provided  in 
the  field  ((Value  of  the  error». 

The  results  of  the  work  are  created  by  pushing  the  button  ((Form  files».  The  exit  from  the  DirTask 
program  is  made  by  pushing  the  button  ((Finish)). 

DielScan  software  environment  serve  from  solution  of  inversed  problems  of  electromagnetic  prob- 
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ing  of  dielectric  structures.  The  Program  has  a  standard  graphic  user  interface  for  programs,  working  under 
Windows  OS.  Upon  program  start,  there  appears  a  interface  window  on  the  monitor,  shown  in  Fig.3.  This 
window  contains  several  different  components.  In  the  upper  part,  immediately  under  the  name  of  the  win¬ 
dow  there  is  a  line  with  a  list  of  the  main  menu  of  the  program.  Below  the  line  of  main  menu  there  is  the 
main  toolbar  as  a  set  of  buttons,  which  duplicate  the  most  frequent  commands  of  the  main  menu. 


^  Solution  of  inverse  problem  of  probing  layered  dielectrics  -  DielScan  lfl?l[X| 


File  Edit  View  Parameters  Work  Information 

a#  e  ►  ■  n  m  k?  f 

STATISTICS  OF  OPTIMIZATION  PROCESS 


Iteration#  ,•  52 

FRecord  0  0444940975840023 
Nwrotuslsste  173  5 

Numb,  of  local  minimizations  5 


Time  lapse 
00:00:00 


B  WXRSHSStaOSlSS  Unchecked 
□  HXRfCSSaMtel  Goo_d_  _ 


flporpec  99 

;TBTyc  Final  sMsctealisa 


Fig.3.  General  outlook  of  the  graphic  interface  of  the  DielScan  program 

In  the  lest  part  of  the  interface  window  there  is  schematic  image  of  the  studied  structure,  the  upper 
part  of  which  shows  the  name  of  file,  in  which  the  input  data  for  the  DielScan  program  are  written.  The 
right  part  of  the  interface  window  shows  some  parameters  that  describe  the  optimization  process.  After 
finishing  program  session,  the  results  of  calculation  are  output  in  the  lower  part  of  the  window. 

Commands  of  menu  «File».  Commands  of  this  menu  allow  opening,  saving  and  printing  out  input 
data  for  the  DielScan  program.  Name  of  the  file,  containing  the  input  data  is  provided  by  the  command 
«Load  task».  In  order  save  changes  in  input  data,  made  with  the  help  of  other  commands,  the  commands 
«Save  task»  and  «Save  As»  are  used.  The  first  command  saves  changes  in  the  same  file,  the  second  -  in  the 
file  with  another  name.  The  command  «Save  As»  opens  a  standard  dialogue  window  for  choice  of  the  file¬ 
name. 

Commands  of  menu  «Editor».  Commands  from  this  menu  are  temporary  not  available. 

Commands  of  menu  «View».  Commands  of  this  group  are  simple  switch  keys.  They  switch  on  or 
off  display  regime  of  the  toolbar  and  status  bar  of  the  program. 

Commands  of  menu  «Parameters».  Commands  of  this  menu  allow  editing  values  of  complex  re¬ 
flection  factor.  Upon  performance  of  the  command  «Frequency  parameters))  there  opens  additional  inter¬ 
face  window,  separated  into  tow  parts  (Fig.  4).  In  the  upper  part  of  the  window  there  is  a  spreadsheet  with 
real  and  imaginary  approximate  values  of  complex  reflection  factor  and  estimations  of  their  precision  for 


2/9/2010 


STCU  PROJECT  P-281  - FINAL  REPORT  FF  PAGE  41 

the  chosen  range  of  frequency  changes.  Content  of  table  cells  (approximate  values  of  reflection  factor  pa¬ 
rameters)  of  can  be  edited  with  the  use  of  standard  Windows  means.  With  the  help  of  the  button  «Add  fre- 
quency»  and  «Exclude  frequency))  the  range  of  frequencies  can  be  modified.  In  the  lower  part  of  the  win¬ 
dow  there  is  graphic  image  of  the  dependence  of  reflection  factor  parameters  on  the  frequency. 


Measurements  data  for  reflection  factor.  rO 


Re 

Im 

/n/n 

W 

Approximate  value 

Mm 

Max 

Approximate  value 

Min 

Mew 

0 

0  5000000 

0  471 434426048781 1 

•0  4962467642618754 

-0  4489851676655063 

0  2267324640053341 

02159356800050800 

0  2386657515845621 

1 

08700000 

-0  430342961954197! 

•0  4529925915307342 

-0  4098504399563786 

0  2393551530842992 

0  2279572886517135 

02519527927203150 

2 

1  2400000 

0  373495523762298! 

-0  3931531829076830 

-0  3557100226307609 

0  2421609273777471 

0  2306294546454733 

0  2549062393449568 

3 

1  6100000 

-0  3564137424638851 

-0  3751723604883001 

-0  3394416594894143 

0  2486729541204581 

0  2368313848766272 

0  2617610043373249 

4 

1  9800000 

0 3023970355644921 

-0  3183126690152547 

-0  2879971767280876 

0  241514226120661; 

02300135486863440 

0  2542255011796434 

5 

23500000 

0  292851834550605! 

-0  3082650890006372 

-02789065090958146 

0  2401653226577521 

02287288787216693 

0  2528056027976345 

6 

27200000 

0  244012584063476 

-0  2568553516457650 

-0  2323929372033111 

0  2141998701699001 

0  2039998763522865 

0  2254735475472640 

7 

3  0900000 

-0  21 00667478201 92- 

n  i  niiinnmmnnft- 

-0  2211228924423078 

-0  2000635693525642 

0  2083154451619781 

01983956620590268 

0  2192794159599770 

Fig.4.  Outline  of  auxiliary  graphical  interface  of  the  DielScan  program 


Commands  of  menu  «Work»  serve  to  control  with  work  of  the  DielScan  program. 


3.  Conclusions 


In  spite  of  a  large  number  of  modem  global  optimization  algorithms  [5-7]  and  their  software  im¬ 
plementations  [13],  the  problem  of  development  of  universal  software  program  for  solution  of  inverse 
problems  in  very  topical.  This  is  stipulated  by  the  fact  that  they  do  not  account  specific  character  of  input 
data,  which  is  its  approximateness.  The  proposed  interval  method  of  branches  and  limits  gives  an  possibil¬ 
ity  to  account  measurement  errors,  does  not  need  additional  numerical  parameters,  which  determination  is 
complicated  (e.g.  Liepschiz  constant).  Numerical  experiments  have  shown  the  efficiency  of  DielScan 
software,  which  realizes  this  method  for  solution  of  inverse  problem  of  dielectric  structures  probing.  At  the 
same  time  we  should  note  the  necessity  of  further  updates  and  studies  aimed  on  widening  a  spectrum  of 
solvable  problems  and  integration  of  this  software  with  existing  measuring  complexes. 
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V.  ADAPTATION  OF  THE  THEORY  OF  INVERSE  PROBLEM  SOLUTION  TO  ESTABLISH 
THE  PARAMETERS  OF  THE  MULTI-LAYERED  STRUCTURE  TAKING  THE  ANTENNA’S 
GEOMETRY  INTO  ACCOUNT 

1.  Improvement  of  processing  of  measured  data  at  automated  measuring  bench 
of  microwave  diagnostics 

Range  of  application  and  resolution  capability  of  microwave  method  directly  depend  on  precision 
and  width  of  frequency  range  of  measurement  of  complex  reflection  coefficient.  Therefore  we  have  made 
the  work  of  improvement  of  ultrahigh-frequency  path  of  the  measuring  test  bench  and  algorithms  for  proc¬ 
essing  of  measured  data. 

Thanking  to  the  decrease  of  parasitic  re-reflecting  between  components  of  UHF  path  (Fig.  1)  we 
have  succeeded  to  widen  the  measuring  frequency  range  up  to  band  of  18  GHz  (58. ..76  GHz)  without  sub¬ 
stantial  loss  in  precision  of  the  CRC  measurement. 

For  a  new  frequency  range  we  have  built  calibration  graphs  by  extinction  and  frequency.  The 
bridge  circuit  has  been  calibrated  and  parameters  of  doubled  T-bridge  have  been  determined.  Frequency 
dependence  of  elements  of  dissipation  matrix  is  shown  in  Fig.  2.  Fig.  3  presents  the  root-mean- square  de¬ 
viations  of  measured  signal  from  the  modeled  under  reference  terminations. 

Let  consider  the  antenna  in  the  equivalent  circuit  of  the  test  bench  (Fig.  1)  as  a  face-to-face  linear 
quadrupole.  If  account  all  elements  of  antenna  dissipation  matrix,  the  relation  between  the  measured  CRC 
V(f)  and  CRC  of  flat-layered  enviromnent  Vs(f)  (Fig.  1)  can  be  written  as 


-4 1 - ' - ' - ' - i - 

55  60  65  f  QHz  70  75 

Fig.  2.  Frequency  dependence  of  components  of  dissipation  matrix  of  the  bridge  scheme. 
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Fig.  3.  The  root-mean-square  deviation  of  measured  and  modeled  signals  of  the  bridge  scheme  on  the  fre¬ 
quency. 


V  =  s,  + — or  v  = - 1 - , 

1-S3F5  S1S3-S2-S3V 

where  Si,.., S3  are  complex  factors,  which  can  be  determined  experimentally.  The  Si  is  the  parameter  of 
antenna  CRC  during  its  emission  into  free  space,  S2  is  a  coefficient  of  transmission  of  the  probing  signal 
into  environment,  S3  describes  the  re-reflecting  of  wave  between  antenna  and  environment. 

In  order  to  determine  parameters  of  antenna,  the  plane  metal  screen  has  been  used  as  a  reference 
termination  in  free  space  Vset(f  x).  The  position  of  the  screen  was  changed  from  zero  plane  x  =  0  to  20  mm 
with  a  step  of  Ax=0.2  mm.  On  the  basis  of  measured  CRCs  of  reference  termination  of  the  „antenna-metal 
screen”  system  the  frequency  dependence  of  antenna  parameters  has  been  determined  (Fig.  4).  Fig.  5. 
shows  the  measured  CRC  of  the  screen  VSei(J)  for  different  positions  x  with  accounting  of  antenna.  Geo¬ 
metric  parameters  of  antenna:  distance  from  antenna  aperture  to  zero  plane  85  mm,  aperture  25x25  mm, 
antenna  length  from  nose  to  aperture  is  65  mm. 


4 


.4 _ i _ 1 _ 1 _ 1 _ 1 _ i _ 1 _ 1 _ 1 _ 

58  60  62  64  66  68  70  72  74  76 

f,  GHz 

Fig.  4.  Frequency  parameters  dependence  of  the  model  of  horn  antenna  Si  ...3(f) 


2/9/2010 


STCU 


PROJECT  P-281  - FINAL  REPORT 


PAGE  45 


4 


4 _ i _ i _ i _ i _ i _ i _ i _ i _ i _ 

53  60  62  64  66  68  70  72  74  76 

f,  GHz 


Fig.  5.  Measured  CRC  of  metal  screen  recalculated  for  free  space  Vset(f),  drift  of  screen  by  Ax=0.8  mm 


Let  us  consider  the  multi-layer  dielectric  structure  and  its  reflection  coefficient.  Let  assume  that 
there  is  a  plane  electromagnetic  wave  propagating  in  the  free  space  and  it  falls  perpendicular  onto  the  infi¬ 
nite  (M-l)-layered  environment  on  a  semi-infinite  substrate  with  flat  interfaces,  which  CRC  depends  on 
the  thickness  of  layers  and  electro-physical  parameters  of  materials  (Fig.  6). 

where  hm,  sm,  tgSm  are  thickness,  relative  dielectric  permeability  and  dielectric  dissipation  of  the  mth  layer, 
respectively.  The  CRC  of  flat-layered  environment  Vs(f)  is  determined  in  the  following  way.  The  expres¬ 
sion  for  the  complex  reflection  coefficient  can  be  written  as  [5]: 


YS  = 


Wp-Y, 

Wo+Y, 


Flere  Yl  is  determined,  depending  on  the  number  of  layers,  by  recurrent  relation: 


v  Ym+i(e*V<lm+Y)  +  Wm(exvqm-l)  v  _u/ 

1m  -  VV m  77  “  ;  XM  ~  yr  M  ’ 

H/m(expqm  + 1)  +  Ym+1  (exp  qm  -1) 

where  Wm  =^]e'm/\i0  ;  qm  =  2  j  •  2nf  ■  hm  (1  -  sin2  0) ;  e'm  =  s0sm(l  -j  ■  tg8m)  is  complex 

dielectric  permeability  of  the  /77th  layer;  e0,  Po  is  absolute  dielectric  and  magnetic  permeability  of  vacuum. 


1 

1 


lx 

Fig.  6.  Model  of  multi-layer  dielectric  structure. 
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Experiments  on  measuring  of  CRC  of  metal  screen  at  different  distances  from  the  zero  plane  have 
been  conducted.  Fig.  7  shows  the  dependence  of  the  module  of  measured  CRC  of  metal  screen  \  VS\  at  the 
distance  to  zero  plane  x.  As  results  of  energy  losses  due  to  dissipation  beyond  limits  of  the  antenna  aper¬ 
ture  projection  the  efficient  reflection  coefficient,  determined  by  the  plane  wave  model,  decreases  with  a 
distance.  At  the  solution  of  inverse  problem  this  would  lead  to  incorrect  interpretation  of  parameters  of 
layered  structure,  particularly  in  the  case  of  substantial  length. 


Fig.  7.  Dependence  of  reflection  coefficient  model  on  the  distance  from  the  screen: 

°  experiment;  —  approximation  |F,|~e_a\ 

It  can  be  seen  from  Fig.  7,  the  change  in  reflection  coefficient,  and  hence  the  wave  amplitude  can 
be  approximately  described  by  exponential  dependence  on  the  coordinate  by  the  same  law  as  for  amplitude 
of  plane  wave  in  dielectric  with  losses.  Thus,  losses  of  energy  for  dissipation  can  be  accounted  in  frames  of 
plane  wave  by  introduction  of  additional  dimensionless  parameter  y,  which  does  not  depend  on  coordinates 
and  parameters  of  environment.  For  dielectric  under  condition  of  tgS,„«l,  y«l  the  dependence  of  pro¬ 
gressive  wave  amplitude  on  coordinate  with  account  of  dissipation  can  be  written  as  [1]: 


/  i  _ 

\Am I ~ exp  271/ •  Ven,e0Mo  ' (tg8m  +  Y) •  * 


where  S0 ,  J-Iq  are  absolute  dielectric  and  magnetic  permeability  of  vacuum.  Then  the  complex  dielectric 
permeability  of  the  mth  layer  with  account  of  losses  for  lateral  dissipation: 

e'm  =eo£m(l-./-(tgSw+Y)). 

2.  Experimental  study  of  inverse  problem  for  diagnostics  of  flat-layered  dielectric 

The  inverse  problem  (IP)  is  to  determine  the  vector  of  the  sought-for  parameters  of  flat-layered 
environment  p*  by  the  measured  frequency  dependence  of  the  CRC  [2,3]  and  can  be  solved  by  the  method 
of  functional  minimization 


p*  =  inf  F(p),  F(p)  =  fy* {f.)  -  V‘h  (p ,/;■) 2 , 

PeP  m1 

where  Nf  is  the  number  of  frequency  responses. 

The  CRC  of  fluoroplastic  was  determined  and  the  inverse  problem  on  parameters  of  flat-layered 
structure  p=(h  1,^2,82^62),  where  h\  is  a  distance  to  the  plate,  hi,  £2,  tgd2  -  are  thickness  and  electro¬ 
physical  parameters  of  fluoroplastic.  Fig.  8  shows  measured  and  calculated  by  the  plane  wave  model  CRC 
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for  fluoroplastic  plate.  The  y  parameter  has  been  estimated  experimentally.  If  to  account  the  y  dissipation 
parameter,  parameters  of  dielectric  are  much  closer  to  reference  data  that  those  determined  by  the  plane 
wave  model  (see  Table  1). 

Fig.  9  shows  the  particular  case  of  IP,  namely  the  objective  function  for  determination  of  geomet¬ 
ric  parameters  of  structure  -  position  and  thickness  of  the  plane  F(h\,  h2).  Fig.  10  presents  the  objective 
function  from  dielectric  permeability  and  thickness  of  the  plate  F{h2,  s2).  It  can  be  seen  from  the  figure, 
that  the  objective  function  is  multiextremal,  therefore  methods  of  global  minimization  should  be  used  in 
order  to  solve  it. 

Table  1.  Estimation  of  dielectric  (fluoroplastic)  parameters 


Environment  parame¬ 
ter 

Model  of  plane  wave 

Model  accounting  the 
dissipation 

Additional  data 

y 

- 

(4. 7±0.1)xl0~3 

- 

h2,  mm 

7.78 

8.193  ±0.005 

- 

S2 

2.110 

2.068  ±0.01 

2.07 

tgS2 

5.18x10- 3 

(0.48  ±0.1) xl0~3 

0.4x1  or3 

60  62  64  66  68  70  72  74  76 

f,  GHz 

Fig.  8.  CRC  of  fluoroplastic  plate: 

—  experiment; ""  the  calculated  CRC  for  plane  wave. 


Fig.  9.  Objective  function  of  the  IP  for  determination  of  location  and  thickness  of  the  plate  F(h  ?,  \\2), 

p*=(-0.00249;  0.00819). 
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3.5 

3 

2.5 

f  2 

1.5 


4  5  6  7  8  9  10  11  12 

h2<  m  x  10'3 

Fig.  10.  Objective  function  of  the  inverse  problem  for  determination  of  parameters  of  dielectric  plate 

F(h 2, £2),  global  minimum  p*= (0.0081 9,  2.07). 

3.  Estimation  of  the  influence  of  different  factors  on  results  of  microwave  diagnostics 

Experimental  studies  of  the  influence  of  interfering  factors  on  result  of  diagnostics  of  layered 
structures  have  been  conducted.  One  of  such  factors  is  the  deviation  of  antenna  axis  from  the  normal  to  the 
surface.  Fig.  11.  presents  the  CRC  of  metal  screen  in  the  case  of  1.8°  axis  deviation. 

Fig.  12  shows  the  range  average  reflection  coefficient  and  dispersion  in  relation  to  model  signal 
depending  on  the  angle  between  the  normal  to  the  surface  and  the  antenna  axis. 

It  can  be  seen  from  figures,  the  change  of  angle  between  normal  and  antenna  axis  causes  the  de¬ 
crease  of  the  power  of  the  sensed  signal,  at  angles  below  1 0  the  influence  is  insignificant,  and  then  it  grows 
substantially.  Flowever  such  a  deviation  is  detected  by  sight  and  can  be  removed  without  additional  meas¬ 
urements. 


Fig.  11.  CRC  of  the  screen  in  the  case  of  deviation  of  antenna  axis  by  1.8°  from  the  normal. 
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Fig.  12a.  RC  and  dispersion  for  the  metal  screen  Fig.  12b.  Dispersion  for  dielectric  plate  in  the 
in  the  case  of  tilt  distortion  case  of  tilt  distortion. 


Let  consider  now  the  influence  of  edges  of  the  structure  on  measurements  results.  Fig.  13  presents 
changes  of  the  CRC  during  the  drift  of  the  screen  relatively  to  the  aperture  projection,  Fig.  14  shows  the 
averaged  module  of  the  RC  of  the  screen  long  the  range  during  its  drift. 


1 1 


0.3 


VA/vVv/^ UAa'i"^ 


’CL  ,'vL/ ■  A. 


I kdufr  1  /y' 


58 


60  62  64  66 


70  72  74  76 


Fig.  13.  Reflection  coefficients  from  metal  screen  (a)  and  dielectric  plane  (b),  moving  over  the  antenna:  red 
-  screen  on  the  edge  of  aperture  projection,  blue  -  screed  half-shields  the  aperture  projection,  green  - 

screen  completely  shields  the  aperture  projection. 


Fig.  14.  Dependence  of  RC  module  on  the  position  of  screed  related  to  the  axis  of  antenna  aperture. 
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It  can  be  seen  from  figures,  the  signal  substantially  changes,  when  projection  of  aperture  crosses 
the  edge  of  screen,  during  approaching  of  antenna  to  the  edge  of  the  structure  the  influence  of  edge  can  be 
neglected.  Conclusions  are  related  only  to  the  far-field  region  (d>A/4A),  for  the  near-field  region  addi¬ 
tional  studies  are  needed. 

4.  Modernization  of  the  portable  microwave  NDT  device 

As  was  noted  before,  the  range  of  application  and  resolution  capacity  of  microwave  method  di¬ 
rectly  depends  on  the  width  of  frequency  range  for  measurement  of  complex  reflection  coefficient.  That  is 
why  we  analyzed  possibilities  to  widen  frequency  range  of  the  HK  device.  During  the  preliminary  stage  of 
the  project  we  made  a  decision  to  change  technical  design  of  the  controlled  generator  of  the  3-mm  wave¬ 
length  range.  New  design  foresees  the  use  of  wide-band  generator  controlled  by  voltage  on  the  Gunn  diode 
in  the  8-mm  range  and  frequency  multiplier  for  transition  into  the  3-mm  range.  This  approach  would  pro¬ 
vide  much  wider  (at  least  5  times)  frequency  readjustment  band,  than  the  available  solid  state  oscillators  in 
this  frequency  range. 

The  model  of  the  generator  has  been  assembled  and  tested.  Fig.  1 5  illustrates  appearance  of  the 
model  of  generator  with  frequency  readjustment. 


Fig.  15.  Model  of  the  unit  for  generation  of  3-mm  range. 


The  structural  scheme  of  the  generator  is  presented  in  Fig.  16.  Here  the  voltage-controlled  genera¬ 
tor  (VCG)  is  the  waveguide  structure  with  the  7. 2x3.4  mm  cross-section  of  the  waveguide.  In  the  middle 
of  waveguide  chamber,  perpendicular  to  the  axis  of  the  waveguide  there  are  Gunn  diode  and  varactor  di¬ 
ode.  Both  diodes  are  in  cases.  From  one  side  of  the  chamber  there  is  the  short  circuitor,  from  the  output 
side  there  is  an  iris.  At  the  exit  from  generator  in  order  to  decrease  the  influence  of  the  outside  load  and 
provide  stability  of  regime  of  soft  agitation  of  generator  the  ferrite  isolator  (FI1)  is  installed. 

Current  of  Gunn  diode  in  working  regime  is  ~1.2  A,  feeding  voltage  is  -3.6  V.  Voltage,  supplied 
to  the  varactor  diode,  changes  in  the  range  +0.5  V  ...  20  V.  In  such  a  case  frequency  of  generator  changes 
in  the  range  of  31. ...33  GHz  under  output  power  of  20. ..35  mW. 

Since  the  changes  in  environment  temperatures  change  the  current  of  Gunn  diode,  geometric  di¬ 
mension  of  waveguide  chamber  and  parameters  of  varactor  diode,  under  constant  controlling  voltage  the 
frequency  of  the  output  signal  will  be  changed  on  varactor  diode.  In  Table  2  and  Fig.  17  we  provide  de¬ 
pendencies  of  signal  frequency  from  generator  at  different  environment  temperatures  (it  was  measured  by 
thermal  sensor  on  the  case  of  generator).  In  can  be  seen  from  figures,  the  changes  of  frequency 
Af~30  MHz  at  AT=50°C  (Uvar  =  10  V).  It  is  also  seen  that  the  decrease  of  Uvarthe  temperature  instability  of 
generator  frequency  increases. 

After  the  ferrite  isolator  on  generator  waveguide  attenuator  (AT)  with  fixed  attenuation,  after  that 
follows  attenuator-modulator.  The  structure  made  of  p-and-n  diodes  installed  inside  the  waveguide  section 
is  used  as  controlling  element  of  the  modulator.  It  is  used  for  modulation  of  generator  signal  (modulation 
depth  is  of  -3  dB  and  installed  during  adjustment  of  detectors  unit). 
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The  signal  comes  to  the  power  amplifier  (YM).  The  fixed  attenuator  (AT)  is  used  for  setting  the 
level  of  output  power  of  driving  oscillator  at  the  entrance  of  amplified.  The  optimal  power  level  at  the  en¬ 
trance  of  amplifier  is  in  the  range  1.5. ..2. 5  mW.  The  amplifier  uses  the  MHC  CHA3093a  transistor  mi¬ 
croassembly  (K=  18... 20  dB,  P_i  =  20...22  dBm,  power  supply  ~  350  mW,  supply  voltage  +3.5  V).  With 
the  help  of  the  second  attenuator  (AT)  the  level  of  the  output  power  at  the  entrance  of  frequency  multiplier 
(FM),  Pin=  100... 130  mW  is  provided  (see  Table  2). 

Frequency  multiplier  "x3"  (FM)  is  the  waveguide  structure  of  overlapping  cross”.  Along  the  axis 
of  waveguides  intersection  there  is  a  compound  coaxial  concordance  structure,  which  serves  as  a  filter  as 
well.  From  one  side  of  the  coaxial  structure  there  is  a  diode  with  Schottky  barrier  of  honeycomb  structure. 
The  contact  with  the  diode  is  provided  with  the  help  of  tungsten  needle.  Agreement  of  multiplier  at  en¬ 
trance  and  exit  (wavelength  channel  of  2.4x1. 2  mm2)  is  performed  with  the  help  of  short-circuitor  pistons. 

As  a  results  of  tuning  of  the  FM  under  conditions  of  the  maximum  bandwidth  we  obtained  levels 
of  output  power  Pout>  1  mW,  the  output  signal  frequencies  changing  in  the  range  of  92. 4. ..98. 5  GFlz  (Ta¬ 
ble  3,  Fig.  18). 


Fig.  16.  Overall  and  functional  scheme  of  the  voltage-controlled  generator  block: 

VCG  -  generator  on  the  Gunn  diode;  FI  -  ferrite  isolator;  AT  -  fixed  attenuator;  M  -  attenuator- 
modulator;  PA  -  power  amplifier;  FI  -  ferrite  circulator  (isolator);  FM  -  frequency  multiplier;  SI  -  power 
source  for  Gunn  oscillator  (+12  V,  1.5  A,  controlling  voltage  Uvar=0.5...20  V);  S2  -  power  source  for 

power  amplifier  (+5  V,  0.35  A). 


Table  2.  Dependence  of  output  frequency  and  Gunn  oscillator  power  on  control  voltage  Uvar  and  tem¬ 
perature _ _ _ _ 


uvar,  V 

o°c 

20°C 

50°C 

P  (PA), 
mW 

f,  GHz 

P,  mW 

f,  GHz 

P,  mW 

f,  GHz 

P,  mW 

0.5 

30.923 

169 

30.795 

162 

30.651 

153 

101 

1 

30.965 

170 

30.848 

165 

30.726 

154 

104 

2 

31.080 

173 

30.955 

168 

30.869 

160 

109 

3 

31.163 

176 

31.124 

172 

31.079 

164 

115 

4 

31.291 

173 

31.255 

171 

31.227 

164 

119 

5 

31.463 

162 

31.408 

162 

31.401 

156 

121 

6 

31.625 

164 

31.597 

160 

31.567 

152 

122 

7 

31.790 

170 

31.758 

165 

31.729 

156 

128 

8 

31.914 

175 

31.891 

170 

31.875 

161 

131 

9 

32.050 

177 

32.037 

172 

32.009 

161 

130 

10 

32.180 

180 

32.164 

175 

32.151 

159 

129 

11 

32.302 

180 

32.275 

174 

32.264 

154 

129 

12 

32.428 

180 

32.401 

173 

32.363 

148 

130 

13 

32.539 

175 

32.518 

167 

32.484 

146 

130 
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14 

32.609 

167 

32.596 

156 

32.574 

134 

129 

15 

32.717 

171 

32.701 

158 

32.676 

132 

129 

16 

32.787 

160 

32.768 

150 

32.758 

125 

125 

17 

32.870 

167 

32.852 

150 

32.835 

123 

126 

18 

32.918 

174 

32.918 

162 

32.897 

140 

132 

19 

32.974 

178 

32.964 

168 

32.948 

150 

134 

20 

33.020 

180 

33.009 

172 

32.993 

158 

137 

Table  3.  Dependence  of  generator  frequency  and  power  on  control  voltage  Uvar. 


77  V 

vary  v 

20  °C 

f.GHz 

P.  mW 

00.5 

92.40 

1.7 

01 

92.55 

1.7 

02 

92.85 

2.3 

03 

93.36 

3.9 

04 

93.75 

4.4 

05 

93.42 

4.2 

06 

94.80 

5.7 

07 

95.27 

5.7 

08 

95.70 

7.1 

09 

96.11 

6.1 

10 

96.49 

5.5 

11 

96.82 

5.5 

12 

97.20 

5.5 

13 

97.55 

5.3 

14 

97.79 

3.9 

15 

98.10 

2.9 

16 

98.30 

1.2 

17 

98.55 

0.56 

18 

98.75 

0.42 

19 

98.90 

0.35 

20 

99.03 

0.31 

Fig.  17.  Output  frequency  (a)  and  power  (b)  of  Gunn  oscillator 
at  different  temperatures  of  environment. 
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Fig.  18.  Output  frequency  (a)  and  power  of  generator  set  (b). 


bj 


5.  Conclusions. 

By  improvements  in  the  UHF  path  on  the  automated  test  bench  of  microwave  NDT  we  succeeded 
in  widening  of  the  CRC  measuring  range  up  to  18  GFlz  (f  =  56. .78  GFlz),  precision  of  frequency  determi¬ 
nation  ±0.02  GFlz,  extinction  coefficient  in  range  0..30  dB  with  a  precision  of  ±1%,  complex  reflection 
coefficient  is  determined  with  the  precision  of  0.01  a.  u. 

The  model  for  probing  of  layered  environment  has  been  built.  It  accounts  a  complete  antenna  dis¬ 
sipation  matrix  as  a  path-through  quadrupole,  which  accounts  the  availability  of  re-reflections  of  proving 
signal  between  antenna  and  environment.  An  additional  parameter  is  introduced  into  the  model  of  interac¬ 
tion  between  the  probing  field  and  flat-layered  structure  which  describes  the  amplitude  decrease  of  the 
proving  field  with  a  distance  because  of  energy  losses  for  radiation  outside  antenna  aperture.  This  is  shown 
to  increase  reliability  of  solution  of  inverse  problem  on  determination  of  parameters  of  single-layer  dielec¬ 
tric. 

The  influence  of  such  interfering  factors  as  the  tilt  of  antenna  axis  and  edge  of  the  structure  during 
measurements  in  the  far-field  zone  on  results  of  diagnostics  has  been  checked  experimentally. 

The  model  for  the  generator  module  of  the  NDT  device  has  been  assembled,  its  frequency  range, 
power  and  temperature  instability  have  been  determined.  Under  maximum  bandwidth  during  changes  of 
output  signal  frequency  in  the  range  92.4. ..98. 5  GFlz  the  levels  of  output  power  Pout>  1  mW  have  been 
obtained. 
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VI.  CONDUCTING  OF  THE  PRECISE  PHYSICAL  EXPERIMENTS 
1.  Experiment  description 

The  sketch  scheme  of  experimental  devise  is  provided  in  Fig.  1.  Automated  test  bench  allows 
measuring  CRF  of  loadings  in  the  58. ..76  GFlz  frequency  range.  The  test  bench  works  as  an  interferometer 
on  the  basis  of  double  T-bridge.  The  unknown  loading  V  (specimen  under  control)  is  joined  to  one  of 
shoulders  of  this  bridge,  the  calibrated  loading  Vc  -  into  the  other.  The  output  signal  of  the  bridge-  "un¬ 
balance  signal"  -  is  measured  and  on  the  basis  of  the  bridge  scheme  model  the  complex  reflection  factor  of 
the  unknown  loading  is  determined. 

Loading  is  the  successively  joint  transmitting-receiving  horn  antenna  and  a  structure,  which  CRF 
Vs  is  to  be  determined.  In  the  test  bench  model  the  antenna  is  considered  as  feed  through  linear  quadruple, 
described  by  a  complex  scattering  matrix. 

The  relation  between  the  measured  V(f)  loading  CRF  and  flat-layered  environment  CRF  Vs(f) 
(Fig.  1)  can  be  written  as 

v  - _ vr _ 

s  S,S,-S2-S3V 

where  Si,.., S3  are  complex  factors,  which  can  be  experimentally  determined.  Parameters  of  bridge  circuit 
and  antenna  have  been  determined  during  previous  stages. 


Vc(f) 

a,  (p 
calibrating 
loading 


Fig.  1.  Schematic  model  of  the  test  bench: 

f  is  oscillation  frequency  generator;  Urefis  detector  of  the  reflected  wave  signal;  Vc  is  the  CRF  of  the  cali¬ 
brated  loading;  K2t  is  dissipation  matrix  of  the  doubled  T-bridge;  V(f)  is  CRF  of  loading  under  measure¬ 
ment;  S(f)  is  antenna;  Vs(f)  is  CRF  of  flat-layered  environment  in  free  space. 

Fig.  2  shows  the  scheme  of  experiment  on  diagnostics  of  exfoliation  in  dielectric  plate.  Exfoliation 
was  modeled  by  a  gap  between  two  dielectric  plates.  Paper  sheets  of  the  0.09±0.01  mm  thickness  were  po¬ 
sitioned  between  plates,  so  the  gap  thickness  was  changed  in  the  0. .  .0.5  mm  range. 

The  plexiglass  was  chosen  and  materials  for  layers  of  model  structure.  This  material  is  cheep  and 
available.  Its  optical  transparency  allows  controlling  the  experiment  process.  The  drawback  of  this  materi¬ 
als  is  a  large  dispersion  of  electrophysical  parameters  between  different  samples  due  to  different  chemical 
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composition.  Therefore  before  modeling  multilayer  structure  parameters  of  each  sample  should  be  deter¬ 
mined. 

vwwwvwvwww 


Fig.  2.  Modeling  of  exfoliation  in  dielectric  layer,  scheme  of  experiment 

Let  consider  the  model  of  interaction  between  the  probing  field  and  the  structure  under  control.  Let 
assume,  that  a  flat  electromagnetic  wave  is  spread  in  free  space  and  falls  perpendicularly  onto  infinite  (M- 

l)-layered  environment  on  semi-infinite  support  with  flat  interfaces,  which  CRF  depends  on  thicknesses  of 
layers  and  electrophysical  parameters  of  materials. 

Vs=Vs(f,  Zll,..,/lM-l)£lv!£M,tg8l,...,tg5M).  (1) 

where  hm,  sm,  tgSm  are  thickness,  relative  dielectric  permeability  and  tangent  of  loss  angle  of  m-th  layer, 
respectively.  The  model  also  accounts  energy  loss  for  scattering  beyond  the  antenna  aperture. 

Let  model  dielectric  plate  with  exfoliation  by  M-l=4-layered  structure,  where  the  layers  are: 

1)  air  gap  between  zero  plane  and  structure,  thickness  of  the  layer  hi  can  be  negative,  air  (s/=1, 
tgS/=0,  fixed  parameters); 

2)  first  dielectric  layer,  thickness  h  2  and  parameters  s2,  tgS2  can  be  unknown  or  fixed,  depending 
on  the  purpose  of  experiment. 

3)  exfoliation,  thickness  h3  is  unknown,  parameters  depend  on  filling  and  correspond  to  either  air 
(£3=1,  tgS3=0,  fixed  parameters)  or  paper  (parameters  should  be  determined); 

4)  the  second  dielectric  layer,  thickness  h4  and  parameters  £4,  tg$4  can  be  unknown  or  fixed,  de¬ 
pending  on  the  purpose  of  experiment. 

The  fifth  semi-infinite  slab  is  air. 

The  inverse  problem  (IP)  for  diagnostics  of  structure  consists  in  determination  of  the  vector  of 
sought  parameters  p*  by  the  measured  frequency  dependence  of  CRF  and  can  be  solved  by  minimization 
of  functional 


P 


=  infF(p),  F(p)  = 

per 


N 


zh 


i= 1 


(2) 


where  Nf  is  the  number  of  frequency  readings.  The  special  software  for  solution  of  the  operatore  equation 
(2)  has  been  developed.  It  includes  effective  algorithm  for  solution  of  direct  problem  (1)  and  minimization 
of  objective  function  F(p)  with  the  use  of  global  minimization  method.  The  program  allows  the  researchers 
to  choose  the  vector  of  sought  parameters  and  fix  known  parameters. 


2/9/2010 


STCU 


PROJECT  P-281  - FINAL  REPORT 


PAGE  56 


2.  Estimation  of  dielectric  plate  parameters 

Before  studying  multilayered  structures  the  plate  material  parameters  should  be  determined. 

The  CRF  of  the  plexiglass  plate  was  measured  and  the  inverse  problem  for  determination  of  flat¬ 
layered  structure  p=(h1,h2,(;2,tgd2),  where  hi  is  the  distance  to  the  plate,  h2,  £2,  tgS2  arc  thickness  and  elec¬ 
trophysical  parameters  of  plexiglass,  has  been  solved.  Fig.  3  shows  the  CRF  of  the  plate,  measured  and 
calculated  by  the  model  of  plane  wave.  Experimentally  measured  CRF  of  the  plate  is  shown  by  solid  line; 
the  model  approximation  resulted  from  the  IP  solution  is  shown  by  dashed  line. 
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Fig.  3.  Evaluation  of  the  plexiglass  plate  parameters: 
h]=— 0.358  mm,  h2=10.561  mm,  s 2=2.574 ,  tgS2-6.31  xlCT3 . 


3.  Experimental  study  on  diagnostics  of  exfoliation  in  dielectric  structures 


In  order  to  study  the  sensitivity  of  radio-wave  method  to  parameters  of  dielectric  structures  we 
have  measured  the  frequency  dependences  of  CRF  of  model  structures.  Dependences  of  CRF  on  thickness 
and  material  of  filling  for  exfoliation  in  the  plate  have  been  studied. 

Fig.  4  illustrates  experimentally  measured  frequency  dependences  of  CRF  of  model  structures  with 
different  gap  thicknesses.  Fig.  5  shows  CRF  for  gaps  of  the  same  thickness  but  with  different  filling  (air 
and  paper).  It  can  be  seen  from  Figures,  that  radio-wave  method  in  far  zone  is  sensitive  to  the  thickness 
and  filling  of  exfoliation  in  dielectric.  The  CRF  of  structure  with  exfoliation  between  dielectric  and  metal 
screen  (Fig.  6)  has  been  studied  as  well. 
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Fig.  4.  Modeling  of  exfoliation  in  dielectric  plate. 

Two  plexiglass  plates  of  10,5  mm  thickness,  air  gap,  filled  by  paper  sheets: 
1  -  without  exfoliation,  2-1  sheet  exfoliation,  3-2  sheet  exfoliation; 
4-3  sheet  exfoliation. 
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Fig.  5.  Modeling  of  exfoliation  in  plate.  Two  plexiglass  plates  10,5  mm,  1  -  exfoliation  of  3 -sheet  thick¬ 
ness  with  air;  2  -  exfoliation  with  2  sheets  of  paper. 
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Fig.  6.  Modeling  of  exfoliation  between  dielectric  and  metal,  plexiglass  plate  10,5  mm,  A1  screen;  green  - 
without  exfoliation,  red  -  exfoliation  of  2  sheets,  black-  exfoliation  of  6  sheets. 


4.  Study  of  the  inverse  problem  for  diagnostics  of  exfoliation  in  dielectric 


Let  consider  the  inverse  problem  for  determination  of  parameters  in  the  structure  with  exfoliation. 
Fig.  7... 9  present  experimental  and  theoretically  approximated  dependences  of  CRF  for  structures  with 
different  exfoliation  depth  (h3).  The  measured  dependence  of  CRF  Vsex(f)  is  shown  by  solid  line,  modeled 
dependence  of  CRF  nymcTHpoM  -  Monc.ibOBany  Vsth(f),  obtained  from  the  solution  of  inverse  problem  is 
shown  by  dash  line.  We  can  see  that  the  proposed  model  for  reflection  factor  describes  the  experimental 
data  rather  good. 
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Fig.  7.  Two  plexiglass  plates  without  a  gap 
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Fig.  8.  Two  plexiglass  plates,  gap  of  2  sheets 
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Fig.  9.  Two  plexiglass  plates,  gap  of  3  sheets 


Several  options  of  the  inverse  diagnostics  problems  have  been  considered: 

1)  determination  of  structure  position,  thickness  and  dielectric  permeability  of  exfoliation:  p=(hi,  h3,  £3); 

2)  determination  of  structure  position,  thickness  of  the  first  layer  and  exfoliation  p=(hi,  h2,  h3); 

3)  determination  of  structure  position,  thickness  of  the  first  layer  and  exfoliation  and  the  third  layer:  p=(hi, 
h2,  h3,  h4). 

Thickness  of  exfoliation  was  changed  from  0  to  6  sheets  of  paper  with  0,09±0,01  mm  thickness. 
Results  of  diagnostics  are  compiled  in  Tables  1..3.  The  average  values  and  dispersion  of  estimated  parame¬ 
ters  of  structures  are  provided  as  well.  The  thickness  of  exfoliation  can  differ  from  the  thickness  of  inset 
since  sheets  can  have  own  curvature.  Therefore  we  consider  as  a  parameter  the  changes  of  the  determined 
gap  thickness  after  adding  one  sheet  Zlh3,  which  corresponds  to  the  thickness  of  one  paper  sheet,  used  for 
reaching  exfoliation  thickness. 
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Table  1.  Solution  of  IP  for  p=(hi,  h3,  S3). 


Exfoliation,  number  of  pa¬ 
per  sheets  (~0,09  mm) 

0 

1 

2 

3... 

Average 

(from 

1-6) 

cr 

hi*,  mm 

-0,5293 

-0,5111 

-0,4935 

-0,4856 

- 

- 

h3*,  mm 

7,7e-10 

0,060 

0,155 

0,245 

- 

- 

zlh3*,  mm 

- 

0,0601 

0,0950 

0,0902 

0,0915 

0,003 

£3* 

- 

0,380 

0,678 

0,807 

0,923 

0,065 

Table  2.  Solution  of  IP  for  p=(hi,  h2,  h3) 


Exfoliation,  number  of  pa¬ 
per  sheets  (~0,09  mm) 

0 

1 

2 

3... 

Average 

(from 

1..6) 

cr 

hi*,  mm 

-0,5245 

-0,5080 

-0,4911 

-0,4852 

- 

- 

h2*,  mm 

10,519 

10,528 

10,542 

10,549 

10,540 

0,007 

h3*,  mm 

0,027 

0,0881 

0,1743 

0,2590 

- 

Ah3*,  mm 

- 

0,0613 

0,0862 

0,0846 

0,0854 

0,0008 

Table  3.  Solution  of  IP  for  p=(ht,  h2,  h3,  h4). 


Exfoliation,  number  of  pa¬ 
per  sheets  (~0,09  mm) 

0 

1 

2 

3... 

Average 

(from 

1..6) 

cr 

hi*,  mm 

-0,5245 

-0,5079 

-0,4911 

-0,4849 

- 

- 

h2*,  mm 

10,664 

10,586 

10,570 

10,564 

10,573 

0,008 

h3*,  mm 

0,03559 

0,0961 

0,1854 

0,2736 

- 

- 

Ah3*,  mm 

- 

0,0604 

0,0893 

0,0883 

0,0888 

0,0005 

h4*,  mm 

10,407 

10,490 

10,510 

10,513 

10,504 

0,009 

Fig.  10  displays  the  dependence  of  CRF  for  the  structure  with  exfoliation  filled  by  paper.  The  in¬ 
verse  problem  on  determination  of  structure  position,  thickness  of  exfoliation  and  its  dielectric  permeabil¬ 
ity  has  been  solved.  It  can  be  seen  that  the  solution  for  the  thickness  of  exfoliation  is  close  to  the  thickness 
of  two  sheets  of  paper. 
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Fig.  10.  Two  plexiglass  plates,  3a3op  2  jihcth,  3anoBHemia  nanip, 
h3*=0,179mm,  £,3*=3,066. 
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Fig.  1 1  illustrates  the  CRF  from  plexiglass  plates  and  metal  screen,  between  which  the  exfoliation 
of  different  width  is  created.  The  position  of  the  plate  and  the  thickness  of  exfoliation  have  been  deter¬ 
mined. 
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Fig.  1 1 .  Plate  made  of  plexiglass,  gap  and  A1  screen: 

1  -  h3*=0,021  mm  (without  exfoliation);  2  -  h3*=0,593  mm  (6  plates). 


5.  Conclusions 

For  different  setups  of  inverse  problem  of  diagnostics  the  estimated  structure  parameters  are 
somewhat  different.  This  is  related  to  errors  in  measurements,  deviation  of  actual  geometric  and  electro¬ 
physical  parameters  from  the  set,  and  by  errors  of  the  model,  which  describes  the  interaction  between  field 
and  environment.  In  most  cases  these  discrepancies  are  within  limits  of  errors  of  conventional  thickness 
measurements,  e.g.  by  micrometer. 

For  the  first  problem  (see  Table  1)  the  estimated  dielectric  permeability  of  exfoliation  83*  differs 
from  physically  possible,  however,  with  the  rise  of  the  gap  it  approaches  to  the  real  one  (s=1  for  air).  Simi¬ 
larly,  the  geometric  thickness  of  the  gap  in  the  case  #1  differs  from  experiments  #2  and  #3.  This  can  be 
explained  by  the  fact  that  the  thickness  and  dielectric  permeability  in  the  direct  problem  are  related  pa¬ 
rameters.  It  has  been  shown  experimentally,  that  the  thinner  is  the  layer,  the  more  difficult  is  to  determine 
simultaneously  its  thickness  and  dielectric  permeability. 

Radio-wave  method  in  far  irradiation  zone  with  multiffequency  probing,  measurement  of  CRF  and 
solution  of  IP  allows  measuring  thicknesses  of  separate  layers  of  structure  and  their  electrophysical  pa¬ 
rameters  with  precision,  sufficient  for  practical  application.  It  is  sensitive  to  changes  in  thickness  of  struc¬ 
ture  of  0.01  mm  range  and  exfoliations  in  dielectric  plate  of  0.1  mm  range  of  thickness. 
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VII.  INVESTIGATION  OF  THE  REAL  DIELECTRIC  STRUCTURES  AND  THE  DEVELOP¬ 
MENT  OF  THE  DIAGNOSTICS  TECHNIQUE. 

1.  Development  of  the  method  for  testing  flat-layered  dielectric  structure  for  apparatus  of 
non-destructive  testing. 


On  the  previous  stages  of  the  project  we  have  considered  the  method  of  NDT  of  layered  dielectric 
structures,  based  on  the  determination  of  complex  reflection  factor  (CRF)  from  the  structure  accounting 
properties  of  measuring  tract,  antenna  and  the  use  of  model  of  interaction  between  probing  wave  and  envi¬ 
ronment.  The  CRF  is  an  important  informative  parameter  for  study  of  the  interaction  between  the  probing 
field  and  the  object  under  control.  Such  an  approach  allows  observe  illustratively  the  connection  between 
parameters  of  structure  and  electromagnetic  field,  estimate  adequacy  of  the  model  and  sensitivity  of  the 
method  to  the  changes  of  environment  parameters.  Measurement  of  the  CRF  is  the  main  mean  for  building 
of  mathematical  model  of  measuring  tract  and  transceiver  antenna. 

Flowever,  for  calculation  of  the  CRF  by  the  proposed  method  one  should  measure  output  signal  of 
the  bridge  scheme  at  several  different  (not  less  than  three)  values  of  the  phase  of  the  reference  signal,  so 
the  calibrated  phase-twister  should  be  readjusted  during  measurements.  This  complicates  the  procedure, 
which  is  not  acceptable  for  the  NDT  device. 

Therefore  we  propose  to  use  to  use  the  output  signal  of  the  bridge  scheme  on  many  frequencies  as 
informative  parameters  and  on  the  basis  of  the  model  of  this  signal  to  determine  parameters  of  dielectric 
structure. 

Let  remind  the  model  of  output  signal,  described  in  the  previous  report.  Fig.  1  presents  the  model 


Fig.  1.  Schematic  model  of  test  bench: 
f  is  sweep  frequency  generator; 

Uref  is  detector  of  the  signal  from  reflected 
wave; 

Vc  is  the  CRF  of  calibrated  loading; 

K2t  is  a  matrix  of  scattering  of  the  double 
T-bridge; 

V(f)  is  the  CRF  of  loading  under  meas¬ 
urement; 

S(f)  is  antenna; 

Vs(f)  is  the  CRF  of  the  flat-layered  envi¬ 
ronment  in  free  space 


Detected  voltage  of  the  incident  and  reflected  waves  are  proportional  to  the  power  of  the  signal  in 
the  UF1F  tract  U  ~  P  ~  |Aj2.  It  has  been  shown  the  unbalance  signal  is  determined  by  the  relation: 

U 


D(v,  vc  (a,  <p),  k0 ,  kx ,  k2 ,  k3 )  = 


£0  4-  k] 


V  -V„ 


(i-k3-Vy(i-k3-Vc)-k2-V-Vc 


(i) 


where  k0(f),..  k3(f)  are  complex  factors,  that  depend  on  the  geometry  of  the  2T-joint  and  frequency.  They 
describe  properties  of  interferometer  and  can  be  determined  on  the  basis  of  data  of  natural  experiment. 

The  relation  between  the  measured  CRF  V(f)  and  CRF  of  the  flat-layered  environment  Vs(f)  can  be 
written  as 


S  V 

V  =  S,+ — or  Vs  = 
1  -S3Vs 


st-v 


s1s3-s2-s3r 


(2) 


where  Si,.., S3  are  complex  factors,  determined  through  elements  of  the  antenna  scattering  matrix.  The  CRF 
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of  flat-layered  environment  depends  on  thickness  of  layers  and  electrophysical  parameters  of  materials,  it 
can  be  written  as 


Parameters  of  bridge  scheme  and  antenna  are  determined  at  the  stage  of  calibration  and  are  known, 
which  has  been  described  in  the  previous  report.  So,  one  can  write  the  generalized  model  of  the  output  sig¬ 
nal  of  the  bridge  scheme 

Dth=Dth(f,  p),  (4) 

where  p  is  vector  of  unknown  (sought)  parameters  of  layered  environment,  f  is  a  frequency. 

In  the  previous  report  we  have  shown  that  on  the  basis  of  the  frequency  dependence  of  the  CRF  of 

the  structure  V™  (/] )  its  parameters  have  been  determined.  We  shall  show  below  that  with  the  use  of  sca¬ 
lar  signal  Dex(f)  one  can  also  determine  parameters  of  the  structure. 

The  inversed  problem  (IP)  for  testing  of  environment  is  to  determine  the  vector  of  sought  parame¬ 
ters  of  flat-layered  environment  p*.  We  can  write  the  target  function  by  the  measured  frequency  depend¬ 
ence  Dex(f)  and  model  of  the  signal  of  bridge  scheme.  The  IP  is  solved  by  the  method  of  functional  mini¬ 
mization 


P*  =  in£ ^d(p),  FD{p)  =  YJ\Dex(f)-Dth(f,pf  ,  (5) 

pep  tr 

where  Nf  is  the  number  of  probing  frequencies. 

Fig.  2  illustrates  the  shape  of  measured  signal  for  free  space  and  dielectric  plate  (plexiglass 
10.5  mm). 


a) 
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Fig.  2.  Output  signal  1 


As  can  be  seen  from  the  figure,  the  signal  gas  large  dynamic  range  and  is  rather  rapidly  changing. 
This  is  related  to  the  fact  that  shoulders  of  the  bridge  are  loaded  to  the  structures  of  different  length.  One 
can  see  also  that  the  signal  is  sensitive  to  the  structure  in  free  space,  however  it  is  difficult  to  analyze  it  in 
such  a  form.  Therefore  further  on  we’ll  consider  difference  signal  Dex(f)=  Dex(f)  -  Dexo(f),  where  Dexo(f) 
is  the  signal  for  free  space  (Fig.  2b). 

Fig.  3  presents  the  dependence  of  the  D(f)  signal  on  parameters  of  layered  structure,  the  solid  line 
shows  the  signal  from  the  structure  without  defects  -  two  plexiglass  plates  without  the  air  gap,  the  dashed 
line  is  a  signal  for  “flawed”  structure,  the  air  gap  between  plates  is  about  0.25  mm.  As  can  be  seen  from  the 
figure,  the  signal  is  sensitive  to  changes  of  structure  parameters. 
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Fig.  3.  Sensitivity  of  the  D  signal  to  the  availability  of  air  gap  between  two  dielectric  plates 
10.5  mm:  ( - )  without  air  gap,  ( - )  gap  0.25  mm. 

Let  consider  now  the  result  of  modeling.  Fig.  4  presents  the  measured  and  modeled  signals  D(f)  for 
dielectric  plate  made  of  plexiglass.  Dependence  is  provided  in  larger  scale  of  frequency. 
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b) 


fom  bridge  scheme  for  free  space  ( — )  and  dielectric  plate  ( - )  (a)  and  difference 

of  these  signals  (b). 
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Fig.  4.  The  D(f)-D0(f)  signal  for  dielectric  plate;  ( - )  experiment, 

( - )  numeric  model. 

In  order  to  compare  results  of  the  solution  of  the  IP  let  consider  the  same  model  structure,  as  they 
were  in  the  previous  chapter,  namely  plexiglass  dielectric  plate  and  two  plates  with  air  gap. 

In  order  to  determine  parameters  of  the  material  (plexiglass)  the  CRF  of  one  plate  has  been  meas¬ 
ured.  Proceeding  from  the  solution  of  IP  we  have  determined  the  following  structure  parameters:  plate  po¬ 
sition  hi*=-0.358  mm,  thickness  h2*= 10.461  mm,  and  electrophysical  parameters  S2*=2,624, 
tg82*=6,71x10  3The  thickness  of  the  plate  was  controlled  as  well  with  the  help  of  micrometer.  In  the 
measurement  point  in  is  h2=  10,46  ±0,005  mm. 

Dielectric  plate  with  an  air  gap  was  described  by  the  model  of  4-layer  environment:  air  layer  be¬ 
tween  zero  plane  and  dielectric  of  hi  thickness;  the  second  and  the  fourth  layers  of  dielectric  with  parame¬ 
ters  h2,  h4,  82=64,  tgS2=tgS4  and  the  third  air  layer  -  gap  h3. 

The  inverse  problem  of  probing  of  air  gap  in  plate  was  solved  in  the  following  statement.  Let  con¬ 
sider  parameters  of  material  (plexiglass)  and  gap  (air)  to  be  preset.  On  the  basis  of  measured  frequency 
dependences  of  CRF  we  were  to  determine  positions  and  thicknesses  of  all  layers  of  structure,  i.e.  the  vec¬ 
tor  of  the  sought  parameters  p*=(hb  h2,  h3,  h4). 

Table  1  presents  results  of  the  solution  of  the  IP  on  determination  of  thickness  of  individual  layers 
for  the  plate  and  for  two  plates  with  an  air  gap.  For  comparison,  we  provide  thicknesses  determined  from 
the  CRF  of  structure  (FVs)  and  by  the  proposed  method  from  the  signal  of  bridge  and  target  function  FD(5). 
It  can  be  seen  from  the  table,  the  IP  can  be  successfully  solved  on  the  basis  of  scalar  signal  of  the  bridge 
D(f).  Deviations  from  solutions  of  IP  by  the  FVs  are  within  measurement  error. 

Table  1.  Comparison  of  IP  results  obtained  by  the  use  of  CRF  Vs(f)  and  by  output  signal  D(f) 


1  plate 

2  plates  with  an  air 
gap 

Number 
of  layer 

h,  mm 

Fvs 

h*,  mm 

Fd 

h*,  mm 

Fvs 

h*,  mm 

Fd 

h*,  mm 

hj*,  mm 

-0.358 

-0.373 

-0.4852 

-0.48516 

h2*,  mm 

10.46 

10.461 

10.458 

10.464 

10.464 

h3*,  mm 

0.2701 

0.27013 

h4*,  mm 

10.41 

10.414 

10.414 

2.  Functional  scheme  of  the  device  for  radio-wave  NDT 

Let  consider  the  functional  scheme  of  the  devise  for  radio-wave  non- destructive  testing  which  is 
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Fig.  5.  Functional  scheme  of  the  device  of  radio-wave  NDT. 


UHF  tract  in  the  generator  block,  bridge  scheme,  detector  block,  etalon  and  calibrated  loadings, 
and  measuring-receiving  antenna.  Generator  serves  for  generation  of  oscillatory  e.m.  wave  in  the 
92...  99  GFlz  frequency  range.  The  work  of  the  generator  block  has  been  reviewed  in  details  in  the  previ¬ 
ous  report. 

Block  of  control  and  digital  processing  of  information  is  created  on  the  basis  of  single-crystal 
microcontroller  MSP430F169.  The  block  includes  as  well  the  input  cascade  for  measurement  of  the  signal 
of  detector  section;  synthesizer  of  controlled  voltage  to  control  the  frequency  of  generator;  interface  for 
connection  with  the  PC. 

The  microcontroller  includes  as  well  the  programmed  ROM  for  storage  of  the  program  codes  and 
working  constants,  RAM,  universal  anisochronous  receiver-transmitter. 

The  scheme  of  digital  control  over  generator  frequency  is  realized  on  the  basis  of  12-digit  digital- 
to-analog  converter,  measurement  of  analog  information  is  performed  by  12-digit  ADC,  and  the  input  cas¬ 
cade  has  as  well  an  amplifier  with  gain  constant  controlled  through  the  digital-to-analog  converter.  The 
appearance  of  the  module  is  shown  in  Fig.  6.  Principal  electric  scheme  of  the  module  is  presented  in  Fig.  7. 


Fig.  6.  Block  of  control  and  digital  processing  on  the  model  of  the  NDT  device. 
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Fig.  7.  Principal  electric  scheme  of  the  block  of  control  and  digital  processing  of  information. 

The  outer  feeding  block  serves  to  produce  feeding  voltage,  necessary  for  the  work  of  generator 
block  and  block  for  control  and  digital  processing. 

During  the  reported  period  the  generator  block,  block  of  control  and  digital  processing  and  detec¬ 
tor  section,  as  well  as  the  first  version  of  software  for  microcontroller  and  programs  for  data  processing  on 
PC  have  been  developed  and  tested.  During  the  next  quarter  we  plan  to  finish  the  elaboration  of  feeding 
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block,  manufacture  calibrated  loadings  and  finalize  the  development  of  software  for  the  device. 

3.  Experimental  study  of  detector  block 

Detector  block  is  a  directed  coupler  on  crossed  wave-conductors,  the  main  channel  of  which  is 
loaded  on  the  non-reflecting  (coordinated)  loading,  in  the  side  channel  there  is  detector  section  based  on 
diode  with  Schottky  barrier.  Directed  coupler  and  coordinated  loading  serve  to  coordinate  output  channel 
of  the  bridge  scheme  with  detector  block  for  avoiding  the  reflected  wave  in  the  bridge  scheme. 

Detecting  diode  is  positioned  in  the  section  of  waveguide  video  bandwidth  with  1.2x2. 4  mm  cross- 
section,  in  the  body  of  detector  section  there  is  as  well  the  amplifier  of  video  signal,  fed  by  +12  V  voltage. 
Detector  section  serves  to  detect  UHF  signals  in  the  90...  120  GHz  range  both  constant  and  amplitude- 
modulated,  frequency  band  of  amplification  of  the  detected  signal  is  0.03... 300  kHz.  Fig.  8  presents  the 
appearance  of  detector  section  and  detector  block  assembled.  Fig.  9  illustrates  the  principal  electric  scheme 
of  detector  section. 


Fig.  8.  The  assembled  detector  block.  Fig.  9.  Principal  scheme  of  the  UHF  detector:  D2 

is  Schottky  diode,  Ampl  is  operational  amplifier 
OP27,  Pin2  is  CP50-751  connector 


Main  technical  characteristics  of  the  detector  section: 

Range  of  working  frequencies,  GHz  78...  100 

Transformation  coefficient,  V/W  >  10  000 

Video  bandwidth,  Hz  1  -  1*103 

Feeding  voltage,  V  12  ...  15 

Weight,  g,  not  more  than  60 

Temperature  of  environment,  °C  from  0  to  +50 

Relative  humidity  at  25  °C,  %  80  max. 

Complex  testing  of  detector  section  have  been  performed.  Fig.  10  shows  it  frequency  and  dynamic 
characteristics.  Dynamic  characteristic  is  monotonous  and  close  to  the  linear  one.  Frequency  dependence 
of  transformation  factor  should  be  taken  into  account  during  processing  of  the  measured  signal. 

Fig.  11  shows  frequency  dependence  of  elements  of  scattering  matrix  of  the  directed  coupler.  It 
can  be  seen  from  the  coefficient  of  transmission  from  input  arm  to  detector  section  S3 1  is  lower  than  - 
15  dB,  and  cross  coupling  is  lower  than  -25  dB  in  the  whole  frequency  range 

Let  now  consider  parameters  of  the  ready-assembled  detector  block.  The  main  requirements  to  de¬ 
tector  block  are  sufficient  sensitivity  of  detector  and  amplifier,  as  well  as  maximum  reflection  factor  to 
damp  down  the  reflected  wave  in  the  bridge  scheme.  Fig.  12  illustrates  frequency  dependence  of  transfor¬ 
mation  coefficient  of  the  detector  block.  It  has  certain  heterogeneity  in  frequency  because  of  interaction  of 
incident  and  reflected  waves  of  the  main  and  branched-off  channels.  In  Fig.  13  we  can  see  the  reflection 
factor  Sn  and  coefficient  of  standing  wave  (VSWR)  of  detector  block.  In  the  whole  frequency  range  it  is 
smaller  than  -20  dB,  and  has  substantial  heterogeneity  in  frequency.  This  can  be  attributed  to  superposi¬ 
tion  of  waves  reflected  from  the  detector  and  from  coordinated  loading  in  the  main  channel  of  directed 
coupler. 
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Fig.  10b.  Dependence  of  output  voltage  on  the 
incident  power  of  UF1F  field  at  95  GFlz  fre¬ 
quency. 


Fig.  11.  Coefficient  of  transmission  and  cross  coupling  of  channels  of  directed  coupler. 


Fig.  12.  Frequency  dependence  of  the  output  Fig.  13.  Frequency  dependence  of  the  reflection  fac- 

voltage  of  detector  block  at  the  input  power  of  tor  and  the  factor  of  standing  wave  of  the  detector 

Pin=100  W.  block. 

The  output  voltage  does  not  have  noticeable  valleys  in  the  working  range  of  frequencies,  and  the 
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signal  reflected  from  detector  block  will  extinct  not  less  than  in  100  times  in  comparison  with  the  incident 
one.  Let  draw  a  conclusion  that  the  detector  block  is  manufactured  at  high  quality  and  can  be  used  to  detect 
the  output  signal  of  the  bridge  scheme. 

4.  Development  of  software  for  controlling  module  of  digital  processing  of  information  and 
programs  for  exchange  with  PC 

Software  for  single-crystal  microcontroller  is  the  background-prompt  system,  which  fulfills  the 
following  tasks: 

1.  Connection  with  PC,  receiving  constants  for  the  work  of  the  system.  Work  regimes  for  detector 
diode,  starting  and  ending  frequency  and  the  scanning  step,  time  of  delay  and  time  of  data  collection  dur¬ 
ing  scanning  are  programmed  from  the  PC.  Constants  can  be  stored  in  the  nonvolatile  random-access 
memory  and  loaded  when  power  supply  is  switched  on. 

2.  Adjustment  of  generator  frequency  according  to  the  measurement  program. 

3.  Measurement  of  signal  from  quadratic  detector  proportional  to  the  power  of  the  UHF  field  in  de¬ 
tector  section,  control  of  variable  amplification  coefficient  of  the  input  cascade,  digital  demodulation  of  the 
measured  signal. 

4.  Data  transmission  to  the  PC  for  further  processing. 

Microcontroller  manufactured  by  Texas  Instruments  Incorporated  MSP430F169  has  embedded  12- 
digit  ADC  (workability  of  200  000  symbols  per  second)  and  digital-to-analog  converter.  Microcontroller 
with  the  help  of  digital-to-analog  converter  carries  out  direct  control  of  the  measurement  scheme,  with  the 
help  of  ADC  conducts  measurements  (gets  signal  from  primary  transformer  and  makes  it  demodulation) 
and  transmits  the  measured  data  to  the  PC.  The  microcontroller  is  connected  to  the  PC  by  the  USB  inter¬ 
face.  In  the  ROM  of  microcontroller  there  is  software  developed  under  the  IAR  Systems.  The  software  for 
the  PC  is  developed  in  the  Delphi  7  shell.  For  the  work  with  the  USB-port  from  the  Delphi  7  shell  the  addi¬ 
tional  TXComm  component  and  the  driver  of  the  virtual  COM-port  have  been  used.  The  TXComm  com¬ 
ponent  provides  high-level  communication  with  the  microcontroller  through  the  virtual  COM-port  thank¬ 
ing  to  built-in  realization  of  protocol  of  data  exchange  through  the  port.  In  particular,  there  are  realized 
procedures  of  opening  and  closing  of  the  port,  writing  and  reading  from  it,  tracing  the  status  of  port,  con¬ 
trolling  the  readiness  signal,  work  with  buffer. 

Software  of  the  PC  contains  the  following  procedures: 

StartCommPort  -  start  of  work  with  microcontroller; 

StopCommPort  -  end  of  work  with  microcontroller; 

SetCommPort(b:byte;D:word)  -  transfer  of  the  package  of  data  to  mictocontroller; 
XCommlData(Sender:  TObject;  const  Received:  Cardinal),  i 
GetData(Received:Word;  Var  wAl,wA2:Word)  -  getting  data  from  mictocontroller. 

Data  exchange  protocol  from  the  PC  to  the  microcontroller  foresees  package  transmission  of  data. 
Each  package  from  the  PC  contains  three  bytes:  the  first  one  is  the  attribute  of  the  command  and  its  num¬ 
ber,  and  the  next  two  are  data  bytes.  To  control  the  measuring  scheme  from  PC  the  following  input  pa¬ 
rameters  are  transferred  to  microcontroller: 

1.  Starting  code  of  digital-to-analog  converter  corresponding  to  the  starting  frequency  of  the  UHF  genera¬ 
tor; 

2.  End  code  of  digital-to-analog  converter  corresponding  to  the  end  frequency  of  the  UHF  generator; 

3.  Change  step  of  the  code  of  digital-to-analog  converter  corresponding  to  the  step  of  changes  of  UHF 
generator  frequency; 

4.  Time  of  measurement  at  each  frequency  point,  its  value  corresponds  to  the  number  of  obtained  observa¬ 
tion  for  demodulation  of  the  signal  from  the  UHF  tract  at  each  frequency  point; 

5.  Amplification  factor  of  the  modulated  signal  before  analog-to-digital  conversion.  Codes  0  to  3  corre¬ 
spond  to  amplification  in  1,  10,  100,  1000  times,  respectively. 

After  receiving  input  parameters  the  microcontroller  starts  working  with  ADC,  digital-to-analog 
converter  and  amplifier.  Codes  corresponding  to  the  frequency  of  the  UHF  signal  are  passed  by  turns  to  the 
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digital-to-analog  converter.  During  measurement  time  for  each  frequency  point  the  ADC  performs  certain 
number  of  transformations  by  groups  of  16.  Microcontroller  carries  out  summing  of  each  group  and  oper¬ 
ates  with  them  as  with  individual  observations  and  calculates  results  of  demodulation  and  initiates  their 
transfer  to  the  PC. 

Measurement  results  are  passed  by  microcontroller  in  the  form  of  a  sequence  of  four  bytes.  Two 
high-order  bits  correspond  to  the  number  of  bytes;  six  low-order  bits  contain  six  bits  of  measurement  re¬ 
sults. 


Algorithm  of  data  exchange  between  PC  and  microcontroller 

1.  Opening  of  port,  in  the  case  if  it’s  impossible,  the  error  message. 

2.  Signal  of  readiness  to  receive  data  (DTR)  with  a  tag  brClear  (unreadiness). 

3.  Cleaning  of  output  buffer  of  the  port. 

4.  Transmission  of  scanning  parameters  to  microcontroller,  they  consist  of  7  packages  of  3-byte  length  (3 
variable  of  the  Byte  type).  The  first  byte  corresponds  to  the  name  of  parameter,  two  the  next  two  bytes 
contain  its  value,  composed  from  16-digit  number  of  the  Word  type. 

The  starting  code  of  the  digital-to-analog  converter,  ending  code  of  digital-to-analog  converter  and 
a  step  of  code  change  of  digital-to-analog  converter,  delay  time  for  collecting  measurement  result  and  am¬ 
plification  factor  are  set  by  16-digit  numbers.  The  available  values  of  amplification  factor  are  1,10,  100, 
and  1000. 

Correspondence  of  the  value  of  the  first  byte  and  the  name  of  the  parameter: 

3  -starting  code  of  digital-to-analog  converter  (two  bytes); 

4  -  end  code  of  digital-to-analog  converter  (two  bytes); 

5  -  changing  step  of  digital-to-analog  converter  code  (two  bytes); 

6  -  time  of  delay  for  storage  of  measurement  result  (two  bytes); 

7  -  amplification  factor  (two  bytes) 

5.  Cleaning  of  input  buffer  of  the  port. 

6.  Delay. 

7.  Sending  the  signal  of  readiness  to  receive  data  (DTR)  with  tag  brSet  (ready). 

8.  Waiting  for  data  receiving. 

9.  Microcontroller  receives  measured  data  from  12-digit  ADC,  results  of  measurements  are  16-digit  num¬ 
bers,  each  of  which  is  the  sum  of  16  observations  at  certain  frequency  point.  To  pass  results  into  the  port, 
4-byte  packages  are  formed  (2  bits  in  each  byte  correspond  to  its  number)  and  subsequently  passed  to  the 
port. 

10.  PC  traces  the  state  of  the  port  and  if  there  are  not  less  than  four  bytes  it  chooses  4  bytes  and  checks 
there  sequence  by  the  control  sum. 

1 1 .  If  the  control  sum  mismatches,  data  are  shifted  until  reaching  the  needed  sequence. 

12.  Upon  matching  of  control  sum  from  received  four  bytes,  a  16-diit  number  is  formed,  which  is  the  result 
of  measurements. 

13.  Measurements  results  are  transmitted  to  the  video  monitor  and  saved  in  a  file. 

14.  After  finishing  of  measurement  the  port  is  closed. 

5.  Conclusions 

A  new  method  for  control  of  parameters  of  layered  dielectric  environment  without  need  of  read¬ 
justment  of  phase-twister  during  measurement  and  calculation  of  environment  CRF  has  been  proposed. 
This  allows  substantially  simplify  the  design  of  the  device,  shorten  the  measurement  time  and  improve 
work  operability  of  the  NDT  system.  The  testing  of  this  method  revealed  that  its  accuracy  and  resolution 
are  not  worse  than  those  in  the  method,  based  on  the  determination  of  CRF. 

The  works  on  modernization  of  the  NDT  device  have  been  continued.  The  detector  block  has  been 
manufactured  and  tested.  The  module  of  control  and  digital  processing  of  measured  information  and  soft¬ 
ware  for  the  device  to  control  the  generator,  measure  and  detect  signal  of  detector  and  transfer  of  informa¬ 
tion  to  the  PC  have  been  elaborated. 


2/9/2010 


STCU  PROJECT  P-281  - FINAL  REPORT  FF  PAGE  72 

VIII.  VERIFICATION  OF  THE  RESULTS  OBTAINED 


1.  Finishing  of  development  of  the  radio-wave  NDT  device  and  software 

During  the  last  quarter  the  design,  calibration  and  testing  of  the  NDT  device  for  flat-layered  dielec¬ 
tric  structures  have  been  finished.  Fig.  la  presents  the  appearabnce  of  the  device  in  the  body  during  run  of 
experiment.  In  the  figure  there  are  also  power  unit,  horn  aerial,  sample  holders,  cables  for  connection  with 
power  unit  and  PC. 

Fig.  lb  presents  UHF  tract  of  the  device  from  the  bridge  circuit  side.  In  the  centre  of  the  photo 
there  is  a  doubled  T-bridge,  which  is  joined  to  the  calibrated  direct  reading  phase  changer  (in  the  left),  and 
horn  aerial  (in  the  right)  and  detector  section  of  the  input  shoulder  of  the  bridge  (at  the  top).  The  doubled 
T-bridge  contains  in  its  body  90°  turns  of  wave-conductors  of  cross-arms,  this  allowed  mo  7jo3bohhjio 
BHBecTH  flanges  of  wave  conductors  for  attaching  reference  loads  onto  the  side  panel  of  the  device  without 
additional  adapters.  The  detector  section  is  connected  to  the  microprocessor  by  coaxial  screened  cables 
with  plug-and-sockets. 

Calibrated  phase  changer  (Fig.  2)  is  designed  as  a  piston,  moving  along  wave-conductor  channel 
under  the  action  of  spring-loaded  micrometer  screw,  this  allowing  control  the  movement  with  0.1  mm  pre¬ 
cision.  Knowing  wavelength  in  the  wave-conductor  one  can  find  the  phase  change  angle  of  the  reflected 
wave. 

The  power  unit  (PU)  for  the  device  has  been  designed  and  manufactured  (Fig.  la).  It  serves  for 
formation  of  supply  voltage  for  generator  and  microprocessor  blocks  and  for  detector  section.  The  power 
unit  is  connected  to  the  device  by  multiple-core  cable  with  a  header.  PU  is  fed  from  mains  system  by 
-220V  AC. 

Fig.  4  presents  the  appearance  of  the  main  user  interface  window  during  measurements.  The  user 
can  set  frequency  range  of  measurement,  number  of  frequency  readings,  and  data  acquisition  at  each  fre¬ 
quency  and  see  real  time  measurement  results  of  the  output  signal.  There  exists  the  possibility  to  save  cali¬ 
bration  curves  and  output  normalized  data  to  the  screen.  Measurements  results  are  stored  in  text  format  for 
further  use  and  treatment. 


b) 


Fig.  1 .  Appearance  of  the  NDT  device  for 
flat-layered  dielectric  materials: 
a)  device  during  the  measurement; 
b)  bridge  circuit  of  the  device; 
c)  appearance  of  the  device  without  upper 
cover  of  the  body. 
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Fig.  2.  Direct  reading  phase  changer,  cross-section  (a)  and  appearance  (b). 
Fig.  3  provides  the  functional  scheme  of  the  device. 


Table  1.  Supply  voltage  of  the  NDT  device 


Supply  voltage,  V, 
max  current,  A 

Purpose 

9  V,  1.5A 

Supply  Gunn  oscillator 

6  V,  1.0A 

Supply  power  amplifier 

12  V,  0.5A 

Supply  detector  section 

24  V,  0.5A 

Supply  module  for  controlling  varactor  diode 
voltage 

f 

generator  unit 


Ug+9v  Uvar  Ua  Umod 

ph 


-220V 
power  unit 


DAC  1 2bi: 


ase  shifter 


Udet  +12V 


to 


ADC  1 2  b  it 


MSP430F1 69 
microprocessor  unit 


1 

PC 


Fig.  3.  Functional  scheme  of  the  radio-wave  NDT  device. 
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Fig.  4.  Work  of  the  device  software. 

2.  Calibration  of  generator  block  by  frequency  and  extinction 

In  order  to  control  the  generator  frequency  the  microprocessor  block  forms  a  code  for  16-digit 
digital-analogue  converter  (DAC),  the  voltage  from  the  DAC  is  passed  to  the  amplifier  with  constant  coef¬ 
ficient  and  then  -  to  varactor  diode.  In  order  to  calibrate  the  device  by  frequency  the  dependence  of  volt¬ 
age  of  varactor  diode  Uvar  on  the  control  code  has  been  measured  with  the  help  of  digital  voltmeter.  Meas¬ 
urement  results  are  presented  in  Fig.  5.  It  can  be  seen  from  the  figure  that  the  voltage  on  the  varactor  start¬ 
ing  from  1.9  V  changes  linearly  and  does  not  fall  below  this  point  due  to  scheme  protection  of  varactor 
diode  against  current  overloading.  Generation  frequency  varies  in  the  93  -  98.5  GFlz  range  nonlinearly 
because  of  dependence  of  generation  frequency  on  controlling  voltage  Uvar.  Further  on  we  shall  consider 
frequency  dependences  of  measured  parameters. 

Dynamic  and  frequency  characteristics  of  generator  block  and  detector  section  have  been  stidied. 
Fig.  6  illustrates  frequency  dependence  of  signal  of  generator,  received  at  the  output  of  detector  section 
after  digital  demodulation.  Detector  section  is  joined  to  the  generator  block.  As  can  be  seen  from  the  fig¬ 
ure,  the  power  if  signal  strongly  depends  on  frequency  and  has  a  maximum  in  the  middle  of  the  range  and 
quickly  decays  in  high  frequency  range.  This  occurs  because  of  complex  frequency  characteristics  of  gen¬ 
erator  block,  in  particular  strong  frequency  dependence  of  transmission  coefficient  is  revealed  by  fre¬ 
quency  multiplier.  As  a  result  the  error  in  measurements  in  upper  frequency  sub-range  will  be  increased 
because  of  low  power  of  measured  signal,  which  leads  to  decrease  of  capacity  of  ADT  and  temperature 
instability  of  measurements.  Frequency  point  at  97.3  GFlz  should  be  omitted  from  consideration  since 
there  characteristics  of  generator  change  very  quickly  and  in  vicinity  of  that  point  there  could  appear  large 
measurements  errors. 

Fig.  7  shows  dependence  of  demodulated  signal  of  detector  section  on  power  of  incident  signal. 
Measurements  were  conducted  in  all  frequencies  of  the  range  with  the  help  of  polarisation  attenuator.  It 
can  be  seen  from  the  figure  that  dynamic  characteristics  of  measuring  tract  are  close  to  the  linear  one. 


2/9/2010 


STCU 


PROJECT  P-281  - FINAL  REPORT 


PAGE  75 


N 

X 

o 


99 

98 

97 

96 

95 

94 

93 


Code  DAC 


Code  DAC 


Fig.  5.  Calibration  of  the  device  by  frequency,  dependence  of  voltage  on  varactor  (a)  and  frequency 

(6)  on  the  DAT  code. 


10 


10 


Pout,  r.u. 


Fig.  6.  Dependence  of  output  signal  of  detector  on  Fig.  7.  Dependence  of  output  signal  of  detector  on 
frequency.  signal  power. 


.Further  on  the  measured  output  signal  of  bridge  circuit  we  shall  note  as  D  and  present  it  in  relative 
units  or  decibels  in  dependence  on  frequency. 


3.  Calibration  of  the  bridge  circuit 

In  order  to  calibrate  bridge  circuit  parameters  of  2T-bridge  should  be  defined.  The  scheme  of  the 
model  is  presented  in  Fig.  8.  According  to  this  model  the  detected  output  signal  of  the  bridge  circuit  (un¬ 
balance  signal)  is  determined  by  form: 
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Vc(f) 

a,  cp 
calibrating 
loading 


Fig.  8.  Schematic  model  of  test  bench: 
f  -  generator  of  oscillating  frequency; 
Uref-  detector  of  signal  of  reflected  wave; 
Vc  -  CRF  of  calibrated  load; 

K2t  -  dispersion  matrix  of  the  double  T- 
bridge; 

V(f)  -  CRF  of  load  under  measurement; 
S(f)  -  aerial; 

Vs(f)  -  CRF  of  flat-layered  environment 
in  free  space 


d(v,V',K,~A)= 


+  kx  ■ 


V-V 


(\-kyV)-(\-k,-Vc)-k2-V  -Vc 


(1) 


where  k0(f),..  k3(f)  are  complex  coefficient  dependent  on  geometry  of  2T-bridge  and  frequency.  They  de¬ 
scribe  properties  of  interferometer  and  can  be  determined  on  the  basis  of  natural  experiment  results. 

Calibration  of  measuring  scheme  includes  measurement  of  parameters  of  bridge  circuit  k0(f),..  k3(f) 
for  all  frequencies  of  the  range.  To  do  so,  the  calibrated  phase  changes  are  joined  to  cross-arms  of  the 
bridge  and  the  dependence  of  signal  D  on  the  phase  of  calibrated  load  for  all  frequencies  is  measured. 
Value  of  the  phase  depende  of  the  choice  of  reference  plane  -  zero  position  of  the  scale  of  phase  changer. 
The  scale  position  ssh  =  0  mm  has  been  chosen  as  zero  plane  BuyiiKy.  Distance  between  neighbouring  posi¬ 
tion  of  phase  changer  was  ssh  =  0.2  mm.  Phase  of  the  CRF  of  reference  load  sh  was  determined  through 
the  length  of  the  wave  in  wave-conductor  a  by  formula 

(P,/,=27t^+  o,  [rad].  (2) 


where  sSh  is  the  position  of  phase  changer  scale;  o  is  phase  incursion  corresponding  to  zero  plane,  this 
parameter  is  dependent  on  frequency  and  is  determined  together  with  parameters  of  bridge  circuit. 

Coefficients  of  mathematic  model  of  the  bridge  are  determined  by  the  following  objective  function 


Fi\k0,..,k^)=  -D[yet,Vc(ai,(pi\k0,...,kiY 


i=\ 


(3) 


where  D;  are  measured  values  of  unbalance  signal  of  the  set  of  calibrated  loads  Vc(  ,  )  and  known  refer¬ 
ence  load  Vet;  N  is  the  number  of  fixed  position  of  the  phase  of  calibrated  load.  Optimum  values  of  model 
coefficients  are  determined  by  the  minimum  of  the  objective  function: 


(4,...,4)*=infFi(4,...,4) 


(4) 


Fig.  9  shows  experimental  dependences  of  unbalance  signal  for  different  position  of  the  calibrated 
phase  changer.  Number  of  the  curve  corresponds  to  the  piston  shift  by  0.2  mm.  These  dependences  were 
used  for  identification  of  parameters  of  bridge  circuit.  Fig.  10  provides  frequency  dependences  for  com¬ 
plex  parameters  of  the  bridge  circuit  k0(f),..  k3(f)  and  for  zero  phase. 
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Fig.  9.  Frequency  dependence  of  unbalance  signal  for  different  positions  of  phase  changer. 
Number  of  curve  corresponds  to  the  piston  shift  by  0.2  mm. 


Fig.  10.  Frequency  dependences  of  complex  coefficients  of  the  bridge  circuit  model  kj  and  position  of  zero 

phase  o- 

CRF  of  unknown  load  is  determined  on  the  basis  of  frequency  dependences  of  signal  D  at  differ¬ 
ent  positions  of  calibrated  load  Vc : 

V *  =  inf  F  (V) 

2  ’wher e  F2(V)=YlDi-D(v,V,(a,<pl)X,--,k,))  ■  (5) 

i 

Fig.  11  shows  CRF  of  reference  load  under  different  positions  of  piston,  determined  by  formula  (2) 
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with  the  use  of  obtained  values  of  model  parameters.  It  can  be  seen  from  the  figure  that  module  of  CRF  is 
close  to  1,  phase  changes  linearly  with  frequency  and  corresponds  to  theoretical  phase  incursion  by  for¬ 
mula  (2). 


1.2  - 


_  0.8- 
> 

—  0.6  - 
0.4  - 
0.2  - 

0 - 1 - 1 - 1 - i - 1 - 1 

93  94  95  96  97  98  99 

3 


2 

■o 

co 

F  1 
> 


CO 

-1 


2_i_ i_ i_ i_ i_ i 

93  94  95  96  97  98  99 

f,  GHz 

Fig.  11.  Measured  CRF  of  reference  direct  reading  phase  changer.  Number  of  the  curve  corresponds  to  the 

piston  shift  by  0.8  mm. 


4.  Measurement  of  CRF  of  loads  in  free  space.  Determination  of  parameters  for  the  model 
of  emitter-receiver  aerial 


Let  now  address  to  measurements  of  structures  in  free  space.  We  have  joined  a  pyramidal  horn  ae¬ 
rial  to  the  lateral  shoulder  of  the  bridge  through  a  wave-conducting  adapter  with  90°  turn  (see  Fig.  la).  Ae¬ 
rial  parameters  are  the  following:  length  from  orifice  to  aperture  is  42  mm,  aperture  25x30  mm.  samples 
for  measurements  were  placed  on  the  special  table. 

The  CRF  of  aerial  has  been  measured  when  it  was  loaded  to  free  space  and  metal  screen,  placed  at 
different  distances  from  aperture,  (see  Fig.  12). 

Relation  between  measured  CRF  V(f)  and  CRF  of  flat-layered  environment  Vs(f)  can  be  written  as 


V  =  S]  + 


$2  Ks 

1  -S3Vs 


or  Vc  = 


Si  -V 


S^-Si-SiV 


(6) 


where  Si,.., S3  are  complex  coefficients,  determined  through  elements  of  aerial  dispersion  matrix  and  cam 
be  determined  experimentally. 
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Fig.  12.  CRF  of  aerial  with  free  space  and  a  the  screen.  Number  of  the  screen  corresponds  to  a  shift  of 

metal  sheet  by  0.2  mm  along  the  aerial  axis. 


In  order  to  determine  parameters  of  the  aerial  a  flat  metal  screen  has  been  used  as  a  reference  load 
in  free  space  Vset(f,  x).  Its  position  was  changed  in  the  range  x=0...20mm  from  zero  plane  with  a 
x=0,2  mm  step.  Distance  from  aerial  aperture  to  zero  plane  was  35  mm. 

It  has  been  assumed  for  metal  screen  that 

|Vsetth|  =  1,  arg(Vsetth)  =  2-27ix/A,0  +  ,  (7) 


where  x  is  position  of  the  screen  in  relation  to  zero  plane,  0  is  length  of  electromagnetic  wave  in  vacuum. 

Complex  coefficients  Si(f),..,S3(f)  for  each  frequency  can  be  determined  using  the  experimentally 
measured  values  of  Vet(f,  x)  and  theoretical  CRF  values  of  reference  load  Vsetth(f,  x). 


(S^S^SJ  =mfF3{Sl,S2,S3) 


N<i 


Nx 

F,(sl,s„s,)='£\va(x,)-v(s„s2,s„vK,a,(xl) 


i=0 


(8) 


where  x,=i  x.  Figs.  13  and  14  present  frequency  dependences  for  Si(f),..,S3(f)  coefficients  of  the  aerial 
model  and  CRF  of  reference  load  recalculated  for  free  space  according  to  the  Vs(f)  model  of  aerial  (6). 

Let  study  sensitivity  of  measured  signal  to  parameters  of  layered  structure..  Let  carry  out  meas¬ 
urement  of  signal  D  for  two  dielectric  plates  made  of  plexiglass  with  10.5  mm  thickness,  between  which 
there  was  formed  a  slit  of  0  to  0.6  mm  thickness.  Such  a  structure  models  foliation  in  dielectric  material. 
Fig.  15  shows  dependence  of  difference  signal  D  -  Dfs  (where  Dfs  is  a  signal  for  free  space)  on  frequency 
for  different  values  of  slit  between  plates.  It  can  be  seen  from  the  figure  that  signal  is  sensitive  to  changes 
in  thickness  of  foliation. 
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Fig.  13.  Frequency  dependence  of  the  model  of  horn  aerial. 
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Fig.  14.  Measured  CRF  of  metal  screen,  0.2  mm  shift. 
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Fig.  15.  Frequency  dependence  of  unbalance  signal  for  different  thickness  of  foliation  in  dielectric  plate. 

Similarly  the  experiment  for  dielectric  plate  with  metal  support  has  been  carried  out.  Slit  between 
dielectric  and  metal  plate  models  foliation  between  dielectric  layer  and  metal  base.  It  can  be  seen  from  Fig. 
16  that  the  measured  signal  is  sensitive  to  the  thickness  of  foliation  and  it  can  be  used  for  diagnostics  of 
similar  structures. 


f,  GHz 

Fig.  16.  Frequency  dependence  of  unbalance  signal  for  different  thicknesses  of  foliation  between  dielectric 

plate  and  metal  substrate. 
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5.  Technical  specifications  of  the  radio-wave  NDT  device  for  flat-layered  dielectric  structures 


Frequency  range  of  measurement,  GFlz 

93.2....  98.5 

Power  of  generator  module,  mW 

8 

Frequency  of  modulations,  Flz 

1000 

Unevenness  of  generator  along  the  range,  dB 

12 

Transformation  coefficient  of  detector,  V/W 

>  10  000 

Band  of  video  transmission,  Flz 

1  -  1*103 

Max.  number  of  frequency  readings 

65000 

Error  of  frequency  setting,  % 

<0.02 

Dynamic  range  of  measurement  of  signal,  dB 

0...  -40 

Time  for  settling  the  measurement  regime,  min 

20 

Output  interface 

USB 

Dimensions  (without  aerial  and  PU),  mm 

150x170x90 

Weight  (without  aerial  and  PU),  kg 

1.5 

Power  supply 

-220 V,  50W 

Environment  temperature,  °C 

from  +0  to  +50 

Relative  dampness  at  25°C,  % 

80  max. 

6.  Conclusions 

The  designing  of  the  radio-wave  NDT  device  for  dielectric  structures  has  been  finished.  The  cali¬ 
bration  of  the  device  by  frequency  and  extinction  has  been  performed.  The  parameters  of  the  models  of 
bridge  circuit  and  aerial  for  measuring  the  CRF  of  structures.  The  models  are  shown  to  describe  adequately 
the  probing  process  of  tested  structures.  The  dependence  of  measuring  signal  on  foliation  parameters  of 
dielectric  structures. 

Frequency  range,  sensitivity,  measurement  precision,  models  of  measuring  system  and  interaction 
of  probing  field  with  dielectric  flat-layered  environment,  as  well  as  methods  for  solution  of  inverse  prob¬ 
lems  would  allow  solving  a  wide  range  of  practical  problems  of  non-destructive  testing. 
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